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ABSTRACT 


Estimates  of  the  uncertainties  attached  to  full-scale  predictions  of  submarine 
resistance  and  propulsion  based  on  standard  submarine  tow-tank  model  tests  are 
obtained  by  means  of  a  global  uncertainty  analysis.  The  analysis  takes  into  account 
all  the  component  uncertainties,  including  the  uncertainties  associated  with  the 
prediction  procedure  and  the  measurements  performed  both  at  model  scale  and  at 
full  scale,  which  influence  the  overall  uncertainty  of  full-scale  predictions. 

ADMINISTRATIVE  INFORMATION 

This  investigation  was  sponsored  by  PMS450,  Work  Request  WR10170/ AA,  Task 
Area  F1947001,  Work  Unit  No.  1-5080-097. 

INTRODUCTION 

Estimates  of  the  uncertainties  attached  to  full-scale  predictions  of  submarine 
resistance  and  propulsion  based  on  standard  submarine  tow-tank  model  tests  are 
obtained  in  this  study  by  means  of  a  global  uncertainty  analysis.  The  analysis  takes 
into  account  the  uncertainties  associated  with  the  prediction  procedure  and  the 
uncertainties  of  measurements  performed  both  at  model  scale  and  at  full  scale. 

Thus,  the  uncertainty  analysis  developed  in  the  study  takes  into  account  all  the 
component  uncertainties  which  influence  the  overall  uncertainty  of  full-scale 
predictions. 

The  prediction  procedure,  summarized  in  Appendix  A,  entails  both  resistance 
and  propulsion  tests.  Five  primary  model-scale  variables  are  measured  in  these 
model-scale  tests.  These  measured  primary  model-scale  variables  are  the  carriage 
speed  and  the  drag,  which  are  measured  in  both  resistance  and  propulsion  tests,  and 
the  propeller  rpm,  thrust,  and  torque  measured  in  propulsion  tests. 
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The  five  measured  primary  model-scale  variables  are  used  to  determine 
several  "transformed"  model-scale  variables  via  analytical  relations.  These 
relations  are  given  in  Appendix  A.  The  transformed  model-scale  variables  include 
the  residuary-drag  coefficient,  determined  in  resistance  tests,  and  four 
nondimensional  variables  obtained  via  propulsion  tests  :  the  total-drag  coefficient, 
the  advance  ratio,  the  thrust-deduction  factor,  and  the  propulsive  efficiency. 

The  curves  representing  the  advance  ratio,  the  thrust-deduction  factor,  and 
the  propulsive  efficiency  as  functions  of  the  total-drag  coefficient  are  fundamental 
elements  of  the  model-scale  to  full-scale  extrapolation.  The  relations  used  in  this 
extrapolation  are  given  in  Appendix  A.  The  extrapolation  procedure  is  usually 
implemented  for  a  specified  full-scale  speed  or  for  a  specified  full-scale  shaft 
horsepower.  Both  cases  are  examined  in  the  global  uncertainty  analysis  developed 
in  the  study.  Three  other  cases,  in  which  full-scale  predictions  are  obtained  for  a 
specified  value  of  the  propeller  rpm,  thrust,  or  torque,  are  also  examined  for 
completeness. 

The  uncertainty  analysis,  based  on  classical  expressions  for  the  errors  [1]  and 
elementary  differential  calculus,  is  expounded  in  Appendix  B.  The  Fortran-code 
implementation  of  the  expressions  for  the  uncertainties  obtained  in  Appendix  B  is 
given  in  Appendix  C.  Example  input  and  output  files  associated  with  the  Fortran- 
code  are  also  included  in  Appendix  C.  The  global  uncertainty  analysis  developed  in 
Appendices  B  and  C  provides  a  practical  tool  for  estimating  the  uncertainties  of  full- 
scale  predictions  in  terms  of  component  uncertainties  attached  to  model-scale  and 
full-scale  measurements. 

Full-scale  theoretical  predictions  are  ultimately  compared  to  values  measured 
in  full-scale  trials.  The  observed  differences  between  the  full-scale  measurements 
and  the  full-scale  values  predicted  using  model-scale  tests  are  usually  expressed  in 
the  form  of  a  correlation  allowance  in  the  relation  defining  the  drag  coefficient. 
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The  correlation  allowance  accounts  for  aspects  of  the  full-scale  flow,  such  as 
the  hull  roughness,  that  cannot  be  accounted  for  in  model  tests.  The  correlation 
allowance  also  accounts  for  other  limitations  of  the  model-test  to  full-scale 
prediction  procedure,  notably  errors  that  are  systematically  introduced  into  the 
predictions  due  to  limitations  inherent  to  the  prediction  procedure.  Thus, 
systematic  errors  associated  with  the  characteristics  of  the  tow-tank  and  the 
experimental  set-up  used  in  the  implementation  of  the  prediction  procedure  are 
largely  included  in  the  correlation  allowance,  as  is  attested  by  the  fact  that  different 
correlation  allowances  are  used  for  different  tow-tanks. 

Thus,  the  correlation  allowance  largely  accounts  for  the  systematic  (bias) 
errors  associated  with  the  effects  of  the  tank  bottom  and  walls,  the  influence  of  the 
strut  holding  the  model,  the  strain  gauges,  and  the  electronic  equipment.  Therefore, 
as  long  as  no  significant  changes  are  made  in  the  characteristics  of  the  tow-tank,  the 
strut,  the  strain  gauges,  the  electronic  equipment,  and  the  testing  procedure, 
systematic  errors  attached  to  these  aspects  of  the  prediction  procedure  can  largely  be 
ignored  in  the  uncertainty  analysis  (since  they  are  already  included  in  the 
correlation  allowance  to  a  large  extent,  as  was  noted  previously). 

Errors  stemming  from  residual  waves  in  the  tank  and  differences  between 
full-scale  and  model-scale  depth  of  submergence  can  also  be  ignored  in  the 
uncertainty  analysis  to  a  large  extent,  although  these  errors  are  likely  to  vary  with 
the  design  speed,  and  thus  cannot  be  completely  ignored  in  the  uncertainty  analysis. 
Additional  systematic  errors  due  to  geometrical  imprecisions  of  the  model  clearly 
are  model-dependent,  and  thus  cannot  in  principle  be  ignored  in  the  uncertainty 
analysis. 

In  summary,  it  is  proper  to  ignore  most  systematic  (bias)  errors  in  an 
uncertainty  analysis  of  a  consistent  prediction  procedure  because  these  consistent 
errors  are  largely  included  in  the  correlation  allowance  attached  to  the  prediction 
procedure.  This  general  consideration  and  consideration  of  the  substantial 
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difficulties  involved  in  obtaining  reliable  estimates  of  bias  errors  —  more  precisely, 
of  the  effects  of  the  bias  errors  that  are  not  already  included  in  the  correlation 
allowance  —  suggest  that  a  reasonable  practical  way  of  accounting  for  bias  errors  is  to 
simply  increase  the  precision  (random)  errors  by  means  of  a  multiplicative  factor. 
Specifically,  the  bias  errors  of  the  measured  primary  model-scale  variables  are  taken 
equal  to  the  precision  errors  of  these  variables  in  the  analysis  considered  further  on. 
This  practical  approach  is  justified  by  the  extreme  complexities,  and  uncertainties, 
involved  in  obtaining  realistic  estimates  of  the  influence  of  the  bias  errors  that  are 
not  already  included  in  the  correlation  allowance  upon  full-scale  prediction 
uncertainties. 

The  precision  errors  attached  to  the  measured  primary  model-scale  variables 
are  determined  via  a  statistical  analysis  of  the  repeatability  of  model-scale 
measurements.  This  repeatability  analysis  is  presented  in  Appendix  D. 

Results  of  the  repeatability  analysis  presented  in  Appendix  D  and  of  the  global 
uncertainty  analysis  expounded  in  Appendix  B  are  presented  below  for  several  cases, 
with  the  purpose  of  analyzing  the  contribution  of  the  major  component 
uncertainties  which  influence  the  overall  uncertainty  of  full-scale  predictions. 
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RESULTS  OF  UNCERTAINTY  ANALYSIS 


The  uncertainty  analysis  developed  in  the  study  is  applied  to  a  typical  case, 
which  is  defined  below  (and  in  the  input  file  listed  in  Appendix  C).  The  identifying 
numbers  of  the  model,  the  propeller,  and  the  resistance  (EHP)  and  propulsion  (SHP) 
tests  corresponding  to  the  case  considered  here  are 


model  no. 

propeller  no. 

EHP  exp.  no. 

SHP  exp.  no. 

XXXX 

XXXX 

XXX 

XXX 

The  tank-water  density  and  viscosity  are 


density 

viscosity 

1.937  slug/ft3 

1.084X1 O'5  ftVsec 

The  length  and  the  wetted-surface  area  of  the  model,  and  the  diameter  of  the 
propeller  are 


length 

area 

prop  diameter 

22.697  ft 

138.179  ft2 

0.9986  ft 

The  carriage  speed  and  the  drag  in  the  resistance  (EHP)  tests  are 


speed 

drag 

6.0  knots 

46.66  lbs 

The  carriage  speed,  the  propeller  rpm,  the  total  drag  RT  ,  the  tow  force  AR  ,  and  the 
propeller  thrust  and  torque  in  the  propulsion  tests  are 


speed 

rpm 

drag 

tow  force 

thrust 

torque 

18.0  knots 

923.4 

392.64  lbs 

70.1  lbs 

551.0  lbs 

1501.0  in-lbs 
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The  slopes  of  the  curves  representing  the  advance  ratio,  the  thrust-deduction  factor, 
and  the  propulsive  efficiency  as  functions  of  the  total-drag  coefficient  are 
respectively  equal  to 

advance  ratio  thrust-deduction  factor  propulsive  efficiency 
-0.249  0.067  -0.015 

The  length  and  the  speed  of  the  full-scale  submarine  are  taken  as 


length 

speed 

380  ft 

25  knots 

Finally,  the  viscosity  of  sea-water  is  assumed  equal  to  1.282X10"5  ft2 /sec. 

As  was  already  noted,  results  of  the  repeatability  analysis  presented  in 
Appendix  D  and  of  the  global  uncertainty  analysis  expounded  in  Appendices  B  and 
C  are  presented  for  several  cases  for  the  purpose  of  analyzing  the  contribution  of  the 
major  component  uncertainties  which  influence  the  overall  uncertainty  of  full- 
scale  predictions. 
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Case  Ml :  contribution  of  precision  uncertainties  of  model-scale  measurements 


It  is  instructive  to  begin  by  considering  the  full-scale  prediction-uncertainties 
for  the  case  when  only  the  uncertainties  that  directly  stem  from  model  tests  are 
taken  into  account.  In  this  case,  called  Ml  hereafter,  the  correlation  allowance  and 
full-scale  conditions  (i.e.  the  density  and  the  viscosity  of  sea  water,  the  geometry  of 
the  full-scale  ship  and  propeller,  the  full-scale  values  of  the  speed,  the  propeller 
rpm,  the  thrust,  the  torque,  and  the  shaft  horsepower)  are  presumed  known 
without  uncertainty.  The  density  and  the  viscosity  of  tank  water,  the  length  and  the 
wetted-surface  area  of  the  model,  and  the  propeller  diameter  are  also  presumed 
known  without  uncertainty  in  case  Ml.  Furthermore,  model-scale  uncertainties  are 
taken  equal  to  the  precision  (random)  errors  determined  in  Appendix  D.  Thus,  bias 
errors  attached  to  the  measured  primary  model-scale  variables  are  not  taken  into 
account  in  case  Ml.  Case  Ml  corresponds  to  a  comparison  of  successive  model  tests 
within  a  series  of  consecutive  tests. 

The  analysis  of  the  repeatability  of  measurements  of  the  five  primary  model- 
scale  variables  for  standard  submarine  model  testing  given  in  Appendix  D  shows 
that  relative  precision  uncertainties  in  resistance  (EHP)  tests  are  approximately  equal 
to  0.15%  for  the  carriage  speed  and  1.2%  for  the  drag.  These  uncertainties  are 
indicated  in  the  following  table  : 

Precision  uncertainties  of  model-scale  measurements  in  EHP  tests 


speed 

drag 

0.15% 

1.2% 

Appendix  D  indicates  that  the  relative  precision  uncertainties  in  propulsion  (SHP) 
tests  are  approximately  equal  to  0.15%  for  the  carriage  speed,  0.3%  for  the  propeller 

rpm,  1.2%  for  the  drag  RTand  the  tow  force  AR  ,  and  0.5%  for  both  the  propeller 
thrust  and  torque.  These  uncertainties  are  listed  below  : 


7 


Precision  uncertainties  of  model-scale  measurements  in  SHP  tests 


speed 

rpm 

drag 

tow  force 

thrust 

torque 

0.15% 

0.3% 

1.2% 

1.2% 

0.5% 

0.5% 

The  analysis  and  the  Fortran-code  given  in  Appendices  B  and  C  show  that  the 
model-scale  measurement  uncertainties  defined  above  yield  a  relative  uncertainty 
of  the  residuary-resistance  coefficient  CR  approximately  equal  to  6%  .  The 
prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and  UEHP  for  the 
full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the  previously- 
defined  model-scale  uncertainties  are  listed  in  the  next  table  for  five  cases, 
corresponding  to  rows  number  2  to  6  of  the  table.  These  five  cases  correspond  to 
predictions  for  a  specified  value  of  the  full-scale  speed,  rpm,  thrust,  torque,  or  SHP. 


Full-scale  prediction-uncertainties  for  case  Ml 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

n/a 

0.35% 

2.2% 

2.25% 

2.2% 

1.55% 

rpm 

0.35% 

n/a 

2.3% 

2.4% 

2.4% 

1.8% 

thrust 

1.15% 

1.2% 

n/a 

2.25% 

2.95% 

2.5% 

torque 

1.15% 

1.25% 

2.25% 

n/a 

1.25% 

0.95% 

SHP 

0.75% 

0.85% 

1.95% 

0.85% 

n/a 

1.6% 

The  prediction-uncertainties  given  in  the  second  and  third  rows,  which 
respectively  correspond  to  predictions  for  given  full-scale  values  of  the  speed  or  the 
rpm,  of  the  foregoing  table  are  nearly  identical.  The  prediction-uncertainties  given 
m  the  bottom  row  of  the  table,  corresponding  to  predictions  for  a  given  full-scale 
value  of  the  shaft  horsepower,  are  on  the  whole  somewhat  smaller  than  the 
prediction-uncertainties  listed  in  the  other  rows. 
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The  prediction-uncertainties  listed  in  the  foregoing  table  represent  the 
contribution  of  precision  errors  of  model-scale  measurements  when  all  other 
sources  of  errors  (including  bias  errors  of  model-scale  measurements,  uncertainties 
of  the  density  and  the  viscosity  of  tank  water,  and  model-scale  geometrical 
inaccuracies)  are  ignored.  As  was  already  noted,  this  table  of  prediction- 
uncertainties  corresponds  to  a  comparison  of  successive  model  tests  within  a  series 
of  consecutive  tests.  The  contribution  of  uncertainties  of  the  density  and  the 
viscosity  of  tank  water,  the  model  length  and  wetted-surface  area,  and  the  propeller 
diameter,  are  considered  in  case  M2;  and  the  sensitivity  of  prediction-uncertainties 
to  model-scale  bias  errors  is  considered  in  case  M3. 
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Case  M2  :  contribution  of  uncertainties  of  tank-water  properties  and  model-scale 
geometry 

The  uncertainties  of  the  density  and  the  viscosity  of  tank  water,  the  model 
length  and  wetted-surface  area,  and  the  propeller  diameter  are  taken  as  is  indicated 
below. 


Uncertainties  of  tank-water  properties  and  model-scale  geometry 


density 

viscosity 

length 

area 

prop  diameter 

0.1% 

1.5% 

0.1% 

0.5% 

0.05% 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
uncertainties  of  model-scale  measurements  defined  in  case  Ml  and  the 
uncertainties  of  tank-water  properties  and  model-scale  geometry  now  considered  are 
listed  in  the  next  table  for  five  cases  corresponding  to  predictions  for  a  specified 
value  of  the  full-scale  speed,  rpm,  thrust,  torque,  or  SHP  ,  as  for  case  Ml. 


Full-scale  prediction-uncertainties  for  case  M2 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

n/a 

0.35% 

2.3% 

2.35% 

2.3% 

1.65% 

rpm 

0.35% 

n/a 

2.4% 

2.5% 

2.5% 

1.95% 

thrust 

1.2% 

1.25% 

n/a 

2.25% 

3.0% 

2.55% 

torque 

1.2% 

1.3% 

2.25% 

n/a 

1.3% 

1.0% 

SHP 

0.8% 

0.85% 

1.95% 

0.85% 

n/a 

1.6% 
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These  uncertainties  are  not  significantly  larger  than  those  obtained  for  case 
Ml.  Thus,  the  uncertainties  of  the  density  and  the  viscosity  of  tank  water,  of  the 
length  and  the  wetted-surface  area  of  the  model,  and  of  the  diameter  of  the  propeller 
are  sufficiently  small  that  they  have  insignificant  effect  upon  the  prediction- 
uncertainties. 

Case  M3  :  contribution  of  model-scale  precision  and  bias  uncertainties 

As  in  cases  Ml  and  M2,  only  the  contribution  of  model-scale  uncertainties  are 
considered  in  case  M3.  Thus,  the  correlation  allowance  and  full-scale  conditions  (i.e. 
the  density  and  the  viscosity  of  sea  water,  the  geometry  of  the  full-scale  ship  and 
propeller,  the  full-scale  values  of  the  speed,  the  propeller  rpm,  the  thrust,  the 
torque,  and  the  shaft  horsepower)  are  again  presumed  known  without  uncertainty 
for  the  case  M3  now  considered. 

As  was  already  noted  in  the  introduction,  model-scale  bias  errors  are  taken 
equal  to  the  model-scale  precision  errors  determined  in  Appendix  D  and  listed 

previously  for  cases  Ml  and  M2.  The  total  (precision  +  bias)  model-scale 

1  /2 

uncertainties  which  are  considered  in  case  M3  are  then  equal  to  2  times  the 
uncertainties  considered  in  cases  Ml  and  M2. 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
previously-defined  model-scale  uncertainties  are  listed  in  the  next  table  for  five 
cases  corresponding  to  predictions  for  a  specified  value  of  the  full-scale  speed,  rpm, 

thrust,  torque  or  SHP  ,  as  for  cases  Ml  and  M2.  This  table  indicates  that  the 

1  /2 

prediction-uncertainties  for  case  M3  are  equal  to  2  times  the  prediction- 
uncertainties  for  case  M2,  as  one  expects. 
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Full-scale  prediction-uncertainties  for  case  M3 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

n/a 

0.45% 

3.25% 

3.3% 

3.3% 

2.35% 

rpm 

0.45% 

n/a 

3.35% 

3.55% 

3.55% 

2.75% 

thrust 

1.7% 

1.75% 

n/a 

3.2% 

4.25% 

3.6% 

torque 

1.7% 

1.85% 

3.2% 

n/a 

1.85% 

1.4% 

SHP 

1.1% 

1.2% 

2.8% 

1.2% 

n/a 

2.25% 

The  prediction-uncertainties  for  case  M3  may  be  regarded  as  the  uncertainties 
of  the  full-scale  predictions  obtained  using  standard  submarine  model  testing  in  the 
tow-tank  for  a  specified  full-scale  submarine  and  propeller  geometry  and  specified 
full-scale  conditions.  However,  comparison  of  full-scale  predictions  to 
measurements  in  full-scale  trials  introduces  additional  uncertainties.  These 
additional  uncertainties,  called  full-scale  uncertainties  hereafter,  stem  from  the 
uncertainties  in  the  values  of  the  density  and  the  viscosity  of  sea  water,  the 
geometry  of  the  full-scale  ship  and  propeller,  and  the  values  of  the  full-scale  speed, 
propeller  rpm,  thrust,  torque,  and  shaft  horsepower.  The  contribution  of  these  full- 
scale  uncertainties  to  the  prediction-uncertainties  is  considered  in  case  F. 
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Case  F :  contribution  of  full-scale  uncertainties 

All  model-scale  uncertainties  are  ignored  in  case  F,  which  only  considers  the 
contribution  of  full-scale  uncertainties.  Thus,  all  model-scale  variables,  and  the 
correlation  allowance,  are  presumed  known  without  uncertainty  in  case  F. 

The  relative  uncertainties  of  the  density  and  the  viscosity  of  sea  water,  of  the 
length  and  the  wetted-surface  area  of  the  full-scale  submarine,  and  of  the  propeller 
diameter  are  taken  as  is  indicated  in  the  following  table  : 


Uncertainties  of  full-scale  input  variables 


density 

viscosity 

length 

area 

prop  diameter 

1% 

2% 

0.5% 

i% 

0.1% 

The  uncertainties  of  full-scale  measurements  are  considered  in  Appendix  E  .  The 
total  (precision  +  bias)  uncertainties  of  full-scale  measurements  are  taken  as 


Uncertainties  of  full-scale  measurements 


speed 

rpm 

thrust 

torque 

SHP 

0.6% 

0.4% 

3.0% 

0.9% 

0.9% 

hereafter. 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  associated  with  the 
foregoing  full-scale  uncertainties  are  listed  in  the  next  table  for  five  cases 
corresponding  to  predictions  for  a  specified  value  of  the  full-scale  speed,  rpm,  thrust, 
torque  or  SHP  . 
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Full-scale  prediction-uncertainties  for  case  F 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

0.6% 

0.75% 

3.5% 

2.05% 

2.45% 

2.25% 

rpm 

0.75% 

0.40% 

3.4% 

1.85% 

2.05% 

1.85% 

thrust 

1.8% 

1.75% 

3.0% 

3.15% 

4.7% 

4.6% 

torque 

1.05% 

0.95% 

3.15% 

0.9% 

1.8% 

1.55% 

SHP 

0.85% 

0.7% 

3.1% 

1.2% 

0.9% 

0.9% 

The  prediction-uncertainties  for  case  F,  which  only  considers  the  contribution 
of  full-scale  uncertainties  (with  all  other  sources  of  uncertainties  ignored),  are 
roughly  comparable  (although  somewhat  smaller  on  the  whole)  to  the  uncertainties 
for  case  M3,  which  only  considers  the  contribution  of  model-scale  uncertainties 
(with  all  other  sources  of  uncertainties  ignored). 
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Case  MF  :  contribution  of  model-scale  and  full-scale  uncertainties 


The  contributions  of  both  the  model-scale  uncertainties  and  the  full-scale 
uncertainties  considered  in  cases  M3  and  F,  respectively,  are  now  combined.  The 
prediction-uncertainties  Uspeed,  Urpm,  U thrust,  Utorque,  USHP  and  UEHP  for  the 
full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  for  this  case,  called  MF  hereafter, 
are  listed  in  the  next  table  for  five  cases  corresponding  to  predictions  for  a  specified 
value  of  the  full-scale  speed,  rpm,  thrust,  torque  or  SF1P  . 


Full-scale  prediction-uncertainties  for  case  MF 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

0.6% 

0.85% 

4.8% 

3.9% 

4.1% 

3.25% 

rpm 

0.85% 

0.4% 

4.8% 

4.0% 

4.1% 

3.3% 

thrust 

2.45% 

2.5% 

3.0% 

4.5% 

6.35% 

5.85% 

torque 

2.0% 

2.1% 

4.5% 

0.9% 

2.55% 

2.1% 

SHP 

1.4% 

1.4% 

4.15% 

1.7% 

0.9% 

2.45% 

The  uncertainties  for  case  MF  are  larger  than  the  uncertainties  for  either  case 
M3  or  case  F,  as  one  expects. 
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Case  MFC  :  sensitivity  to  variations  in  the  correlation  allowance 


The  prediction-uncertainties  for  case  MF  are  based  on  the  assumption  that  the 
correlation  allowance  is  known  without  uncertainty.  However,  variations  in  the 
value  of  the  correlation  allowance  occur,  due  to  variations  in  the  full-scale 
submarine  that  are  not  accounted  for  in  the  model  tests  (e.g.  variations  in  the  hull 
roughness)  as  well  as  uncertainties  attached  to  both  model-scale  and  full-scale 
variables.  As  is  noted  in  the  introduction,  bias  errors  systematically  introduced  at 
model  scale  and  full  scale  are  largely,  although  not  fully,  included  in  the  correlation 
allowance. 

The  correlation  allowance  is  taken  equal  to  0.0005  for  the  typical  case 
examined  in  the  present  uncertainty  analysis.  Experience  with  the  SSN  688  class 
submarines  indicates  variations  of  the  correlation  allowance  within  a  fairly  broad 
range.  Inasmuch  as  the  contributions  of  model-scale  and  full-scale  uncertainties  are 
already  included  in  the  full-scale  prediction-uncertainties  obtained  in  case  MF,  a 
20  /o  variation  of  the  correlation  allowance  is  considered  here,  i.e.  variations  of  the 
correlation  allowance  CA  within  the  range 

CA  =  0.0005  +/-0.0001 

are  considered  in  case  MFC  .  Thus,  the  prediction-uncertainties  obtained  when  the 
effect  of  a  20%  variation  in  the  value  of  the  correlation  allowance  is  added  to  the 
model-scale  and  full-scale  uncertainties  considered  in  case  MF  is  examined  in  case 
MFC. 


The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP  for  this  case,  called 
case  MFC,  are  listed  in  the  next  table  for  five  cases  corresponding  to  predictions  for  a 
specified  value  of  the  full-scale  speed,  rpm,  thrust,  torque  or  SHP  . 
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Full-scale  prediction-uncertainties  for  case  MFC 


at  given 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

speed 

0.6% 

0.85% 

6.0% 

5.35% 

5.45% 

4.9% 

rpm 

0.85% 

0.4% 

6.0% 

5.4% 

5.5% 

4.9% 

thrust 

3.1% 

3.1% 

3.0% 

4.5% 

6.6% 

6.15% 

torque 

2.75% 

2.8% 

4.5% 

0.9% 

3.2% 

2.8% 

SHP 

1.85% 

1.85% 

4.35% 

2.1% 

0.9% 

2.45% 

The  uncertainties  for  case  MFC  are  larger  than  the  uncertainties  for  case  MF, 
as  one  expects. 
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CONCLUSION 


In  summary,  a  tool  for  estimating  the  uncertainties  attached  to  full-scale 
predictions  of  submarine  resistance  and  propulsion  based  on  standard  submarine 
tow-tank  model  tests  has  been  developed,  using  a  global  uncertainty  analysis,  and 
applied  to  a  typical  case.  The  analysis  takes  into  account  the  uncertainties  associated 
with  the  prediction  procedure  and  the  uncertainties  of  measurements  performed 
both  at  model  scale  and  at  full  scale.  Thus,  the  uncertainty  analysis  developed  and 
applied  here  takes  into  account  all  the  component  uncertainties  which  influence  the 
overall  uncertainty  of  full-scale  predictions. 

Estimates  of  the  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque, 
USHP  and  UEHP  attached  to  the  full-scale  speed,  rpm,  thrust,  torque,  SHP  and  EHP 
have  been  obtained  for  five  cases,  which  correspond  to  predictions  for  a  specified 
value  of  the  full-scale  speed,  rpm,  thrust,  torque  or  SHP.  The  prediction- 
uncertainties  for  a  specified  value  of  the  full-scale  shaft  horsepower  are  on  the 
whole  somewhat  smaller  than  the  prediction-uncertainties  corresponding  to  a 
specified  value  of  the  full-scale  speed,  rpm,  thrust,  or  torque.  Estimates  of  the 
prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and  UEHP  are 
given  for  six  cases,  called  Ml,  M2,  M3,  F,  MF  and  MFC. 

The  prediction-uncertainties  for  case  Ml  represent  the  contribution  of 
precision  errors  of  model-scale  measurements  when  all  other  sources  of  errors 
(including  bias  errors  of  model-scale  measurements,  uncertainties  of  the  density  and 
the  viscosity  of  tank  water,  and  model-scale  geometrical  inaccuracies)  are  ignored. 
Thus,  bias  errors  attached  to  the  measured  primary  model-scale  variables  are  not 
taken  into  account  in  case  Ml,  which  corresponds  to  successive  model  tests  within  a 
series  of  consecutive  tests. 

The  contribution  of  uncertainties  of  the  density  and  the  viscosity  of  tank 
water,  the  model  length  and  area,  and  the  propeller  diameter,  are  considered  in  case 
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M2.  The  uncertainties  of  full-scale  predictions  for  case  M2  are  not  appreciably  larger 
than  those  obtained  for  case  Ml.  Thus,  the  uncertainties  of  the  density  and  the 
viscosity  of  tank  water,  of  the  length  and  the  area  of  the  model,  and  of  the  diameter 
of  the  propeller  are  sufficiently  small  that  they  have  insignificant  effect  upon  the 
prediction-uncertainties. 

The  sensitivity  of  prediction-uncertainties  to  model-scale  bias  errors  is 
considered  in  case  M3.  Model-scale  bias  errors  are  taken  equal  to  the  model-scale 
precision  (random)  errors  considered  in  cases  Ml  and  M2.  Thus,  the  prediction- 

uncertainties  for  case  M3  are  equal  to  2  times  the  prediction-uncertainties  for  case 
M2  as  one  expects.  The  prediction-uncertainties  for  case  M3  may  be  regarded  as  the 
uncertainties  of  the  full-scale  predictions  obtained  using  standard  submarine  model 
testing  in  the  tow-tank  for  a  specified  full-scale  submarine  and  propeller  geometry 
and  specified  full-scale  conditions. 

Comparison  of  full-scale  predictions  to  measurements  in  full-scale  trials 
introduces  additional  uncertainties.  These  additional  full-scale  uncertainties  stem 
from  the  uncertainties  in  the  values  of  the  density  and  the  viscosity  of  sea  water,  the 
geometry  of  the  full-scale  ship  and  propeller,  and  the  values  of  the  full-scale  speed, 
propeller  rpm,  thrust,  torque,  and  shaft  horsepower.  The  contribution  of  these  full- 
scale  uncertainties  to  the  prediction-uncertainties  is  considered  in  case  F  ,  which 
only  considers  the  contribution  of  full-scale  uncertainties  (with  all  other  sources  of 
uncertainties  ignored).  The  prediction-uncertainties  for  case  F  are  roughly 
comparable  (although  somewhat  smaller  on  the  whole)  to  the  uncertainties  for  case 
M3,  which  only  considers  the  contribution  of  model-scale  uncertainties  (with  all 
other  sources  of  uncertainties  ignored). 

The  contributions  of  both  the  model-scale  uncertainties  and  the  full-scale 
uncertainties  considered  in  cases  M3  and  F  ,  respectively,  are  combined  in  case  MF. 
Thus,  the  uncertainties  for  case  MF  are  larger  than  the  uncertainties  for  either  case 
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M3  or  case  F  .  The  prediction-uncertainties  for  case  MF  are  based  on  the  assumption 
that  the  correlation  allowance  is  known  without  uncertainty. 

However,  variations  in  the  value  of  the  correlation  allowance  occur,  due  to 
variations  in  the  full-scale  submarine  (such  as  variations  in  the  hull  roughness) 
that  are  not  accounted  for  in  the  model  tests,  as  well  as  uncertainties  attached  to 
both  model-scale  and  full-scale  variables.  As  is  noted  in  the  introduction,  bias  errors 
systematically  introduced  by  the  prediction  procedure  and  measurements  at  model 
scale  and  full  scale  are  largely,  although  not  fully,  included  in  the  correlation 
allowance.  Inasmuch  as  the  contributions  of  model-scale  and  full-scale 
uncertainties  are  already  included  in  the  full-scale  prediction-uncertainties 
evaluated  in  case  MF,  the  effect  of  a  20%  variation  in  the  value  of  the  correlation 
allowance  added  to  the  model-scale  and  full-scale  uncertainties  determined  in  case 
MF  is  examined  in  case  MFC  . 


The  cases  Ml,  M2,  M3,  F,  MF  and  MFC  are  summarized  below 
Cases  Ml,  M2,  M3,  F,  MF  and  MFC 


Ml 

only  considers  precision  uncertainties  of  model-scale  measurements 

M2 

adds  uncertainties  of  tank-water  properties  and  model-scale  geometry 

M3 

considers  all  model-scale  precision  and  bias  uncertainties 

only  considers  full-scale  uncertainties 

MF 

considers  all  model-scale  and  full-scale  uncertainties 

MFC 

adds  sensitivity  to  variations  in  correlation  allowance 

The  prediction-uncertainties  Uspeed,  Urpm,  Uthrust,  Utorque,  USHP  and 
UEHP  for  a  specified  value  of  the  full-scale  SHP  are  listed  in  the  following  table  for 
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the  six  cases  Ml,  M2,  M3,  F,  MF  and  MFC.  The  uncertainty  UPC  of  the  propulsive 
efficiency  is  also  given  in  the  table 


Full-scale  prediction-uncertainties  for  a  specified  SHP 


case 

Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

UPC 

Ml 

0.7% 

0.8% 

2.0% 

0.8% 

n/a 

1.6% 

1.6% 

M2 

0.8% 

0.9% 

2.0% 

0.9% 

n/a 

1.6% 

1.6% 

M3 

1.1% 

1.2% 

2.8% 

1.2% 

n/a 

2.3% 

2.3% 

F 

0.8% 

0.7% 

3.1% 

1.2% 

0.9% 

0.9% 

n/a 

MF 

1.4% 

1.4% 

4.2% 

1.7% 

0.9% 

2.4% 

2.3% 

MFC 

1.9% 

1.9% 

4.3% 

2.1% 

0.9% 

2.4% 

2.3% 

In  summary,  it  may  be  concluded  that  the  full-scale  prediction-uncertainties 
for  a  specified  SHP  are  approximately  equal  to 

Summary  of  full-scale  prediction-uncertainties  for  a  specified  SHP 


Uspeed 

Urpm 

Uthrust 

Utorque 

USHP 

UEHP 

UPC 

2% 

2% 

4% 

2% 

1% 

2% 

2% 
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APPENDIX  A  :  PREDICTION  PROCEDURE 
Primary  model-scale  variables 

Primary  model-scale  variables  are  determined  from  measurements.  Five  major  primary 
variables  are  measured  :  the  carriage  speed  V ,  the  drag  R ,  and  the  propeller  rps  n ,  thrust  T 
and  torque  Q . 

Transformed  model-scale  variables 

Transformed  model-scale  variables  are  obtained  from  the  measured  primary  variables  by 
means  of  analytical  relations.  The  major  transformed  model-scale  variables  are  the  total- 
drag  coefficient  Ct  ,  the  friction-drag  coefficient  Cf  ,  the  residuary-drag  coefficient  CR ,  the 
advance  ratio  Jy ,  the  thrust-deduction  factor  1  —  t  and  the  propulsive  efficiency  rju .  The 
relations  defining  these  transformed  variables  are  given  below. 

The  total-drag  coefficient  Cj  of  the  model  in  a  resistance  (EHP)  test  is  given  by  the  relation 

„  Rt 

T  =  Jsvy2  (1) 

where  RT  is  the  drag  of  the  model  (without  propeller) ,  p  is  the  density  of  the  tank  water,  S  is 
the  wetted-surface  area  of  the  model,  and  V  is  the  carriage  speed.  The  friction-drag  coefficient 
Cf  of  the  model  is  evaluated  using  the  ITTC  formula 

Cf  =  0.075/ (ClnRn  —  2)2  (2a) 

where  C  ~  0.4342944819 .  The  Reynolds  number  Rn  is  defined  as 

Rn  =  VL/v  (2b) 

where  L  is  the  length  of  the  model  and  u  is  the  kinematic  viscosity  of  the  tank  water.  The 
residuary-drag  coefficient  CR  is  defined  by  the  relation 


Cr  =  Ct  —  Cp  (3) 

where  Ct  and  Cp  are  given  by  (1)  and  (2) ,  respectively.  The  residuary-drag  coefficient  CR  is 
evaluated  at  low  speed,  so  that  free-surface  effects  may  be  neglected. 

Propulsion  tests  yield  the  ideal  resistance  Ri 

Ri  =  Rt  —  A R  (4) 

where  Rt  is  the  drag  of  the  model  without  propeller  and  AR  is  the  change  in  the  tow  force 
due  to  the  propeller.  The  total-drag  coefficient  Ct  of  the  model  in  a  propulsion  test  is  given 

by 

Ct  =  Jsvy2 

The  advance  ratio  Jy  is  defined  as 

Jv  =  V/(nD)  (6) 
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where  n  and  D  are  the  propeller  rps  and  diameter.  The  thrust-deduction  factor  1  —  t  is  given 

by 

1  -t  =  Ri/T  (7) 

where  T  is  the  propeller  thrust.  The  propulsive  efficiency  r//;  is 


V D  = 


VRj 

2irnQ 


(8) 


where  Q  is  the  propeller  torque.  The  three  curves  representing  the  nondimensional  variables 
Jy ,  1  —t  and  t]d  as  functions  of  Cp  are  used  to  determine  full-scale  predictions. 


Full-scale  predictions 

The  superscript  S  identifies  full-scale  variables.  No  superscript  is  used  for  model-scale 
variables.  At  a  given  ship  speed  F5,  the  total-drag  coefficient  Cp ,  the  propeller  rps  ns , 
thrust  Ts  and  torque  Qs ,  and  the  shaft  horsepower  SHPs  are  determined  using  the  relations 
given  below. 

The  friction-drag  coefficient  Cp  is  determined  using  the  ITTC  formula 

Cp  =  0.075/(Cini?f  —  2)2  (9a) 


where  C  ^  0.4342944819 .  The  Reynolds  number  is  defined  as 

=  VsLs/u$ea  (9b) 

where  Ls  is  the  ship  length  and  vsea  is  the  kinematic  viscosity  of  sea  water.  The  total-drag 
coefficient  of  the  ship  is  evaluated  using  the  relation 

C%  =  Csf  +  Cr  +  Ca  (10) 

Here,  Cp  is  the  friction-drag  coefficient,  Cr  is  the  residuary-drag  coefficient  determined  from 
model  tests,  and  the  correlation  allowance  CA  accounts  for  differences  between  the  actual  drag 
coefficient  C and  the  predicted  drag  coefficient  Cf-  +  Cr  . 

The  propeller  rps  is  obtained  from  the  relation 


n 


s  y!_ 

4ds 


(ii) 


where  Ds  is  the  propeller  diameter.  Furthermore,  the  advance  ratio  Jy  is  determined  from 
the  function  Jy(Cr )  obtained  from  model  tests,  with  Cp  taken  equal  to  the  full-scale  value 
Cj,  predicted  by  (10) . 

The  total  drag  of  the  ship  /?|'  and  the  power  required  to  overcome  i?‘)'  are  given  by 

Rp  =  Cp  psea  Ss(Vs)2/2 


EHPS  =  R^VS  —  Cp  psea  Ss{Vsf/2 


where  psea  is  the  density  of  sea  water  and  Ss  is  the  wetted  surface  of  the  ship. 


(12) 
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The  thrust  Ts  exerted  by  the  propeller  is  evaluated  using  the  relation 


rpj  _ 


-  jrfs  -  e”°ss(vsf/2 


(13) 


where  the  thrust-deduction  factor  1—  ts  is  determined  from  the  function  (1—  t)(Cr)  obtained 
from  model  tests,  with  Cj  taken  equal  to  the  full-scale  value  given  by  (10) . 

The  power  provided  to  the  propeller  is 

SHPs  =  psea  Ss{Vs)3/ 2  (14) 

VdVd 

where  the  propulsive  efficiency  is  determined  from  the  function  t]d(Ct)  obtained  from 
model  tests,  with  Cj-  taken  equal  to  the  full-scale  value  Cj>  given  by  (10) . 

Finally,  the  propeller  torque  Qs  is  defined  by 


SHPs 
2ir  ns 


C$  JyD^_ 
4irr]f) 


psea  SS(yS} 2 


(15) 


where  (14)  and  (11)  were  used. 

The  foregoing  relations  yield  values  of  the  shaft  horsepower  SHP $  corresponding  to  a  range 
of  values  of  the  ship  speed  Vs .  A  plot  of  the  speed  V5  as  a  function  of  the  horsepower  SHPs 
is  then  used  to  determine  the  ship  speed  Vs  corresponding  to  a  prescribed  value  of  the  power 
SHPs . 
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APPENDIX  B  :  GLOBAL  UNCERTAINTY  ANALYSIS 

Uncertainties  of  measured  model-scale  variables 

The  previous  relations,  which  define  model-scale  and  full-scale  variables  in  terms  of  five 
measured  primary  model-scale  variables  (speed  V ,  drag  R ,  and  propeller  rps  n ,  thrust  T  and 
torque  Q ) ,  can  be  used  to  determine  the  uncertainties  of  the  transformed  model-scale  variables 
and  of  the  predicted  full-scale  variables  in  terms  of  the  uncertainties  of  the  measured  primary 
variables.  These  analytical  expressions  for  the  uncertainties  of  the  transformed  model-scale 
variables  and  of  the  full-scale  predictions  are  given  below. 


Uncertainties  of  transformed  model-scale  variables 

Expression  (1)  yields 

dCT  dRr  dp  dS  dV 
Ct  Rt  P  S  V 


(16) 


Using  (2)  we  obtain 

dCF  —2  C  dRg  —2  C 
Cp  C\n.Rn  —  2  Rn  \/0.075 

where  C  =  4C2/0.075  ~  10.0593 .  We  thus  have 

dCF 


r—  fdV  dL  du\  fc—  (dV  dL  du\ 

'Mv  +  T  -  v)  = -Vc^(v  +  T  -  t) 


Expression  (3)  yields 


r— — -  (dV  dL  dv\ 

=  (-  +  T-v) 

dCF 


(17) 


dCR  =  Cp  -pr~  ~  CF 

L'T 


cF 


By  using  (16)  and  (17)  we  obtain 

dCR  =  CT  -j-  §)  +  CFVlOCF[f  -  ^)-(2Cr  — CFV/10CF)  ^ 

Thus,  the  uncertainty  of  the  residuary-drag  coefficient  is  defined  by 

“'©Iff)'*©'] 


(18) 


where  SR  stands  for  SRt  and  V  is  the  carriage  speed  for  the  resistance  test  used  to  determine 
Cr . 

Expression  (4)  yields 

dR{  =  dRr  ~  dAR 

The  absolute  uncertainty  of  Ri  therefore  is  given  by 


(19) 
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Expression  (5)  shows  that  the  relative  uncertainty  of  Ct  is  given  by 


Expression  (6)  defines  the  relative  uncertainty  of  the  advance  ratio  Jy  as 


The  relative  uncertainty  of  the  thrust-deduction  factor  1  -t  is  defined  by  (7)  as 


The  relative  uncertainty  of  the  propulsive  efficiency  r/D  is  defined  by  (8)  as 


In  (20)  and  (22)-(23) ,  6Ri  is  given  by  (19) . 

Uncertainties  of  full-scale  predictions 

Expression  (10)  for  the  total-drag  coefficient  yields 

dC-f  =  dCp  +  <1Cr  +  dC a 
which  may  be  expressed  in  the  form 

dCtp  _  Cp  dCp  dCR  +  dC/A 

~cf  =  c$~cf+  c$ 

Expressions  (9) ,  (2)  and  (17)  show  that  we  have 


CS  V  0CF  yVS  +  LS  vsea  J 

(24) 

We  then  have 

d<^T  —y( pLu  dVS  \ 

cf  l  c  ) 

(25) 

where  T  and  Dqv 

are  defined  as 

r=Jiocf  £f 

(26) 

pLv  _  dCR  +  JCa  dLs  dusea 

~  \JlQCp  CSF  LS  "sea 

(27) 

Expressions  (11) ,  (13) ,  (14)  and  (15)  involve  the  nondimensional  coefficients  Jy  ,  1-t  and 
rjD  .  These  coefficients  are  determined  from  curve  fits  (obtained  from  model  tests)  of  Jy  ,  1-t 
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and  Tjo  as  functions  of  Ct  .  Let  A  stand  for  any  one  of  the  three  coefficients  Jy  ,  1— t  and  t)d  * 
The  difference  in  the  coefficient  A  is  given  by 


dk\K+4^dC* 


The  first  term  represents  the  difference  in  A  at  a  given  value  of  Cs ,  and  the  second  term 
defines  the  difference  in  A  due  to  the  uncertainty  of  Cs .  Let  the  first  term  be  written  as  dA 
for  shortness.  The  relative  difference  in  A  may  then  be  expressed  in  the  form 

dA  dA  Cf  /  L,  <nr‘\ 

A  +  dCT  A  \  c  Vs  j 

where  (25)  was  used.  The  relative  differences  in  the  coefficients  Jy  ,  1 —t  and  ;//)  are  then  given 
by 


dJv  ,  ,/Di,  <*FS\  <*(1-*)  .  T,  (nL,  dr,D  (  Lu  dVs\ 

-jy+”  [D°—vr)  1=r+,T  l°c  "Wj  ^+<T  I®0  (28) 


where  aJ,  af  and  av  are  defined  as 


Here,  expression  (26)  for  the  term  T  was  used. 
It  is  also  useful  to  define  the  notation 


<tv  = 


Cf  drjD 
7 Id  dCT 


(29) 


Dv  = 

d,Vs 

Dn  = 

dns 

_  dTs 

„  dQs 

Dp  = 

dSHPs 

(30a) 

Vs 

ns 

TS 

Dq~  QS 

SHPS 

Dr  = 

dpsea 

psea 

dSs 
+  ss 

Df 

_d(l-t) 

"  1  ~t  V 

=  Dd,  , 

VD 

M 

^  4r 

dDs 
+  DS 

(30b) 

Expressions  (11) ,  (13) , 

(14), 

(15) ,  (25) ,  (28)  and  (30)  yield 

Dn 

=  (l+aJ)Dv- 

aJD j?  -  D f 

(31a) 

Dt  = 

■■(2- 

T  +  (rt)£>v  +  (r- 

-o*)D&  +  Df 

-Dt 

(31b) 

Dq- 

=  (2-r  +  afl- 

■aJ)Dv  +  (T-a 

r’  +  aJ)D^  +  Ds-D 

v  +  Dj 

(31c) 

Dp  = 

(3- 

r  +  cr7')Dv  +  (r- 

-a11)  D(f  +  Ds 

Djj 

(31d) 

where  the  relative  differences  dJy/Jy ,  d(  1—  ts)/(l-ts)  and  drjf) /t]q  have  been  taken  equal 
to  the  corresponding  model-scale  values  dJy/Jy ,  d(  1—  t)/(  1—  t )  and  d-qp/po  ■ 

The  four  relations  (31)  involve  model-scale  variables  and  differences  —  which  occur  via 
expressions  (18) ,  (21) ,  (22) ,  (23)  for  the  terms  dCR,  dJy/Jy  ,  d(l-  t)/(  1—  t) ,  dpo/VD  — 
and  full-scale  variables.  The  full-scale  variables  include  the  relative  differences  dpsea/psea , 
dvseajvsea  ?  dDs j pS ,  dLs  /  Ls  and  dSs/Ss  (which  may  be  determined  independently  and 
thus  may  be  presumed  known  for  the  purpose  of  this  uncertainty  analysis) ,  the  difference  dCU 
(which  may  also  be  regarded  as  a  given  input  for  this  analysis) ,  and  the  five  terms  dVs /Vs, 
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dns /ns ,  dTs /Ts,  dQs /Qs ,  dSHPs / SHPs .  Thus,  four  of  these  five  terms  may  be  determined 
from  any  one  of  them.  Specifically,  the  four  relations  (31) ,  which  define  the  relative  differences 
dns /ns ,  dTs  /Ts,  dQs /Qs  and  dSHPs  /  SHPs  in  terms  of  dVs /Vs,  can  be  expressed  in  the 
four  alternative  forms  given  below. 

One  alternative  form  of  the  four  relations  (31)  is 

( 1  +aJ)Dv  =  DN  +  +  D $  (32a) 

( 1  +  aJ )  Dt  =  ( 2  -  T  +  o* )  (Dn  +  D$ ) 

+{T-ot+2aJ)D%'  +  (1  +  oJ  ){Dsp  -  Dt)  (32b) 

{l+<TJ)DQ  =  (2-r  +  ar’-aJ)DN  +  (r-ar>  +  3<jJ)Dj? 

+  (l+aJ)(Dsp-Dv)  +  (3-T  +  ^)D^  (32c) 

( 1+  aJ )  DP  =  ( 3  —  T  +  </>  ){Dn  +  DDj  )  +  ( T-  +  3 a3 )  D& 

+(l+aJ)(Ds-D.v)  (32d) 

These  relations  define  the  relative  differences  dVs /Vs,  dTs/Ts,  dQs /Qs  and  dSHPs /SHPs 
in  terms  of  dns /ns . 

The  four  relations  (31)  can  similarly  be  expressed  in  the  form 

(2-r +  o*)Dv  =  DT-(r-at)Dt?-Dsp+Dt  (33a) 

( 2  -  T  +  a* )  Dn  =  ( 1+  aJ  )(Dt  -  Ds  +  Dt )  -  ( T-  af+  2 a3 )  D& 

— ( 2  —  T  +  cr* )  Dj  (33b) 

( 2  —  T  +  <r4 )  Dq  =  ( 2  —  T  +  cr7^  —  o"J )  (Dt  +  Dt ) 

+(**-  a”  +  aJ )( 2 Dhcv  +  Dsp  )  +  ( 2  -  T  +  a4  )(D^  -  Dv )  (33c) 

(2-T  +  at)DP  =  ( 3  -  T  +  <tv)(Dt  +  Dt)  -  (T+  2^  -  3cr4 )  D& 

-(l+a7^ -at)Dsp-(2-T  +  at)Dr,  (33d) 

These  relations  define  the  relative  differences  dVs /Vs ,  dns/ns,  dQs/Qs  and  dSHPs /SHPs 
in  terms  of  dTs  /Ts . 

A  third  alternative  form  of  the  four  relations  (31)  is 

{2-r  +  ^-aJ)Dv  =  DQ-(r-^  +  aJ)D^-Dsp+Dri-Df  (34a) 

( 2  —  T  +  a71  -  aJ )  Dn  =  (1+  aJ  ){Dq  -  Ds  +  Dv )  -  ( T-  a71  +  ZaJ )  D% 

-(3-T  +  ar,)DI?  (34b) 
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(34c) 


( 2  -  r  +  a*  -  aJ )  Dt  =  ( 2  -  T  +  a* )  (Dq  +  Dv  -  Dd3  )  -  (a1  -  cr*  +  a' 1 ) ( 2 D%  +  Dsp ) 

-(2-T  +  (7r>-(TJ)Dt 

(2  -  T  +  a*  -  <tj  )  DP  =  (3  -T  +  (t^^Dq  -  D$ )  -  (T-  a*  +  3aJ )  D{? 

-  (1+ aJ )(DS  -  Dv)  (34d) 

These  relations  define  the  relative  differences  dVs/Vs,  dns/ns,  dTs /Ts  and  dSHPs/SHPs 
in  terms  of  dQs /Qs . 

Finally,  the  four  relations  (31)  can  be  expressed  as 

{3  -  r  +  a11)  Dv  =  Dp  -  {r-  a")  Dj?  -  Dsp  +  Dv  (35a) 

(3 -r  +  o»)  DN  =  (1+  aJ)(DP  -  Dsp  +  D„)  -  (IT-  ^  +  3aJ )  D& 

-{3-r  +  av)Df  (35b) 


(3  -  r  +  6* )  Dt  =  (2  -  r  +  a1  )(DP  +  Dv)  +  (T+  2^  -  3<r* )  Dg 

+  (l+^-at)Dsp-(3-T  +  ^)Dt  (35c) 


( 3  -  T  +  o* )  Dq  =  ( 2  —  T  +  a1'  -  aJ )  DP  +  ( T  —  a”  +  3  a' 7 )  Dj? 

+  (l+aJ)(Dsp-Dv)  +  (3-r  +  a^)D ?  (35d) 

These  relations  define  the  relative  differences  dVs /Vs,  dns/ns,  dTs/Ts  and  dQs/Qs  in  terms 
of  dSHPs/SHPs. 

Expressions  (12)  and  (14)  show  that  the  relative  difference 


Dp  = 


dEHPs 

EHPS 


(36) 


can  be  obtained  from  the  foregoing  expressions  for  the  relative  difference  Dp  =  dSHPs  / SHPs . 
Specifically,  expressions  (31d) ,  (32d) ,  (33d)  and  (34d)  respectively  yield 


De  =  {3-Y)Dv  +  TD^ +  Dsp 

( 1+  <7J )  De  =  (3  -  T)(DN  +  DDj  )  +  (T+  3aJ)  D&  +  ( 1+  aJ )  Ds 
(2  -T  +  af)  De  =  (3  -T)(T>t  +  Dt)  -  (T-  3af)  D&  -  (1- a*)  Ds 
(2-T-aJ)DE  =  (3-r)(DQ-D^)-(T+3aJ)D^-(l+aJ)Dsp 
The  relations  EHPS  =  rjE  SHPS ,  (28)  and  (35a)  yield 


(37a) 

(37b) 

(37c) 

(37d) 


(3  -  r  +  o»)De  =  (3  -  r  )(Dp  +  Dv)  +  0*  {3D%  +  Ds )  (37e) 

The  relative  uncertainties  6VS /Vs ,  6ns/ns,  6TS/TS,  6QS/QS  and  6SHPs /SHPs  corre¬ 
sponding  to  the  relative  differences  dVs/Vs,  dns  jns,  dTs/Ts,  dQs/Qs  and  dSHPs/SHPs 
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defined  in  the  foregoing  alternative  relations  are  readily  determined  by  taking  the  square  root 
of  the  sum  of  the  square  of  every  term  in  these  relations.  Thus,  we  define  the  notation 


ubv 


{6CrI±  (6Ca)2  ,  (6LS\2  ,  fSvsea\2 

10(C|)3  +  J +V  J 


UV  = 


( 6SHPs \2 
\SHPS  ) 


UE  = 


( SEHPs  \2 
\EHPS  ) 


(38a) 


(38b) 


(38c) 


Uv  = 


(38d) 


corresponding  to  (27)  and  (30) .  We  also  define  the  relative  uncertainties  attached  to  the 
full-scale  measurements  of  Vs,  ns,  Ts,  Qs,  SHPs  and  EHPS,  i.e. 


ulsm  = 


ufNsm  = 


f6n%m^\ 


TV- 


ufm  = 


(38e) 


u}Psm  = 


(MHPjL)2 

\  SHPS  ) 


fWHPk]2 

EHPS  ) 


(38f) 


where  the  subscript  or  superscript  fsm  means  full-scale  measurement.  In  practice,  the  un- 
certainty  U^m  attached  to  the  effective  horsepower  EHPS  of  the  ship  is  not  defined  because 
measurements  of  EHP^sm  are  not  available.  Thus,  the  term  EHP^sm  in  the  expressions  given 
below  may  be  ignored  in  practice. 


Expressions  (31)  and  (37a)  yield 


UN  =  (l+aJ)2  Ulsm  +  (aJ  )2  USU  +  Uf  +  U}™  (39a) 

UT  =  {2-T  +  ct?U$m  +  {T-ot)2Ubv  +  US  +  Ut  +  Ubsm  (39b) 

UQ  =  ( 2  -  T  +  <7*  -  )2  Ulsm  +  ( T -o"  +  aJ?Uhv  +  U*  +  Uv  +  U?  +  Uq™  (39c) 

Up  =  ( 3  -  r  +  ^  )2  ulsm  +  ( r—  ^  )2  ubu  +  +  uv  +  ufPsm  (39d) 

Ue  =  ( 3  -  r  )2  ulsm  +  r2  Ubv  +  usp  +  UfEsm  (39e) 

Similarly,  (32)  and  (37b)  yield 

( 1+  cr7  )2  UV  =  UfNsm  +  (<rJ  )2  Uhv  +  Uf  +  ( l+aJ  )2  Ulsm  (40a) 


(l+<TJ)2UT  =  (2-r  +  at)2(uftm  +  U?)  +  (r-at+2<TJ)2US1' 

+  (1  +<rJ)HuH  +  Ut  +  U}m)  (40b) 
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{1+  <rJ  )2Uq  =  {2-T  +  a7*  -  aJ  )2U$m  +  {T-  a7'  +  3aJ  )2Uj? 

+  ( 3  —  T  +  cr7*  )2  Uf  +  (\+  crJ  )2  ( Up  +Urj  +  Uq™  )  (40c) 

{l+aJ)2UP  =  {3  -  T  +  a*  )2  {Ufim  +  U? )  +  (T-  a71  +  3aJ  )2  U£u 

+  (l+aJ)2(Uj;  +  Ur)  +  UsPsm)  (40d) 

{l+aJ)2UE  =  {3-r)2{u£m  +  U?)  +  (r+3aJ)2USl/ 

+  (l+aJ)2{U*  +  ugm)  (40e) 

Expressions  (33)  and  (37c)  yield 

( 2  -  r  +  (7*  )2  Uv  =  ufm  +  {r-at)2u£u  +  U$  +  Ut  +  {2-r  +  (Tt)2  ulsm  (41a) 

( 2  -  r  +  a* ): 2  UN  =  ( 1  +  a' 1 )' 2  ( ufim  +  Usp  +  Ut )  +  ( T -  a' 1  +  2  a' 7 ): 2  U£u 

+  (2-r  +  at)2(U§>  +  U$m)  (41b) 

( 2  —  T  +  a*  )2  Uq  =  ( 2  -  T  +  o»  -  aJ  )2  (  Ufm  +  Ut )  +  (<r*  -  a7’  +  a3  )2  ( 4  U&  +  Usp  ) 

+  {2-V  +  atf{Uf  +  Uv  +  Ugm)  (41c) 

( 2  -  r  +  a*  )2  C7p  =  ( 3  -  r  +  cr7'  )2  ( ufm  +  Ut)  +  (F+  2a71  -  3^  )2  u£u 

+  ( 1+  a7*  -  a*  f  U%  +  ( 2  -  T  +  )2  ( Uv  +  Ugm )  (41d) 

( 2  -  r  +  a* )' 2  uE  =  ( 3  -  r ) 2  ( u£sm  +  ut )  +  ( r  -  3  a* )‘ 2  u£v 

+  {\-<jt)2U% +  {2-T  +  at)2Ugm  (41e) 

Expressions  (34)  and  (37d)  yield 

( 2  -  T  +  a71  -  a3  )2  Uv  =  Ugm  +  ( T-  a71  +  aJ  )2  U£u  +  U$  +  Uv  +  Uf 

+  (2-T  +  a71  -aJ)2Ulsm  (42a) 

(2  —  T  +  av  —  aJ  )2  Un  =  (1+  a-7  )2  (  Uq771  +  Up  +  Uv)  +  (T—  av  +  3aJ  )2  U£u 

+  (3  -  r  +  a*  )2uf  +  (2  -  r  +  a71  -  aJ  )2  uftm  (42b) 

( 2  -  T  +  a7'  -  aJ  )2  UT  =  ( 2  -  T  +  a*  )2  ( ugm  +  Un  +  U§> )  +  (a* -  a7*  +  aJ  )2  (4U£V  +  U% ) 

+  {2-T  +  a7t -aJ)2{Ut  +  u£sm)  (42c) 

(2  —  T  +  an  —  crJ  )2  Up  =  {3-r  +  a7>)2{Ugm  +  U?)  +  (r-aT> +  3<rJ)2U£u 

+  ( 1+  a3  )2  ( Usp  +  Uv )  +  ( 2  -  T  +  a71  -  a3  f  u£m  (42d) 
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(2  -  r  -  <jJ)2  uE  =  ( 3  —  r  )2  (  Uq™  +  u?)  +  (t+3<tj)2u&' 

+  (l+aJ)2Uj; +  (2-r-aJ)2UfEsm  (42e) 

Finally,  (35)  and  (37e)  yield 

( 3  -  r  +  a*  )2  Uv  =  Ufpsm  +  ( r -  o*  )2  uSv  +  USP  +  Uv  +  ( 3  -  r  +  a77  )2  ulsm  (43a) 

( 3  -  r  +  a*  )2  UN  =  (1+  crJ  )2  ( UfPsm  +  Usp  +  Uv )  +  (T-  a*  +  ZaJ  )2  U%v 

+  (3-T  +  ari)2(U?  +  Uftm)  (43b) 

( 3  -  r  +  ^  )■ 2  uT  =  ( 2  -  r  +  <7*  )■ 2  ( ufpsm  +  uv )  +  ( r  +  2  ^  -  3  ^  )■ 2  u£u 

+  ( 1+  a77  —  ct*  )2  Up  +  ( 3  —  r  +  a7'  )2  ( +  UEm )  (43c) 

( 3  -  r  +  a77  )'■ 2  Uq  =  ( 2  -  r  +  a77  -  a' 7 ): 2  USPsm  +  ( T -  a77  +  3  a' 7 ): 2  U%v 

+  (l+*J)2(Uj;  +  Ur])  +  (3-r  +  aT>)2(U?  +  Ugm)  (43d) 

(3  -  T  +  a77  )2  UE  =  (3  -  T)2  ( Ufpsm  +  Ur,)  +  (a77)2  (9 u&  +  U%) 

+  (3-r  +  c*)2ugm  (43e) 


The  five  sets  of  alternative  expressions  (39)  through  (43) ,  and  expressions  (38) ,  (29)  and 
(26) ,  define  the  uncertainties  attached  to  the  full-scale  predictions  of  the  ship  speed  Vs,  the 
propeller  rpm  Ns,  thrust  Ts  and  torque  Qs,  the  shaft  horsepower  SHPs  and  the  effective 
horsepower  EHPS  for  five  cases,  in  which  Fs,  Ns,  Ts,  Qs  or  SHPs  are  held  constant  (within 
the  accuracy  of  full-scale  measurements) . 
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APPENDIX  C : 

FORTRAN-CODE  &  INPUT-OUTPUT  FILES 

The  source  file  UA3.f  of  the  program  UA3  ,  which  represents  the  (third 
version)  Fortran  implementation  of  the  uncertainty  analysis  expounded  in 
Appendix  B  ,  is  given  in  Appendix  C.  The  symbols  defined  in  the  analysis  and  used 
in  the  Fortran-code  UA3  are  fairly  consistent. 

An  example  of  the  input  file  UA3.in  required  by  UA3.f ,  and  of  the 
corresponding  output  file  UA3.out  generated  by  UA3.f  ,  is  also  given  in  this 
Appendix.  The  attached  example  input  file  UA3.in  and  output  file  UA3.out 
corresponds  to  the  previously-defined  case  MFC  ,  in  which  the  uncertainties 
attached  to  model-scale  and  full-scale  variables  and  to  the  value  of  the  correlation 
coefficient  are  included. 
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FORTRAN  CODE 


c  General  uncertainty  analysis  of  full-scale  resistance 
c  and  propulsion  predictions  using  tow-tank  model  tests 
c  Francis  Noblesse  (May  97 ;  modified  June  97) 
c 

program  UA3 
c 

character  TTxp*50  ,  date*50  ,  model*50  ,  prop*50  , 

&  EHPxp*50  ,  SHPxp*50  ,  comment*80 
c 

real  ro  ,  nu  ,  Uro  ,  Unu  ,  L ,  S  ,  D  ,  UL  ,  US  ,  UD  ,  VCRknot , 
&  RTCR  ,  URTCR  ,  Vknot ,  RT  ,  DelR  ,  nrpm  ,  T  ,  Qinlb  , 

&  UV ,  URT  ,  UDelR ,  Un ,  UT ,  UQ ,  JVCT ,  tdCT ,  etaCT , 

&  nusea  ,  Uros  ,  Unus  ,  LS  ,  VSknot ,  ULS  ,  USS  ,  UDS  , 

&  UVfsm  ,  UNfsm  ,  UTfsm  ,  UQfsm  ,  UPfsm  ,  CA  ,  dCA  , 

&  VCR ,  V ,  VS ,  n ,  Q ,  CTCR ,  CFCR ,  CR ,  Ri ,  CT ,  CF , 

&  JV ,  td  ,  eta  ,  Uro2  ,  Unu2 ,  UL2  ,  US2  ,  UD2  ,  URTCR2  , 

&  UV2  ,  URi ,  URi2  ,  Un2  ,  UT2  ,  UQ2  ,  UCA  , 

&  UCTCR2  ,  UCTCR ,  UCFCR2 ,  UCFCR ,  dCR2  ,  UCR , 

&  UCT2 ,  UCT ,  UCF2 ,  UCF , 

&  UJV2  ,  Utd2 ,  Ueta2,  UJV ,  Utd  ,  Ueta  , 

&  CFS  ,  CTS  ,  UroS2  ,  UJD2  ,  cofCFS ,  UCLnu2 , 

&  Gamma  ,  sigmaj  ,  sigmat ,  sigmeta  , 

&  UVfsm2  ,  UNfsm2  ,  UTfsm2  ,  UQfsm2  ,  UPfsm2  , 

&  UN2V  ,  UT2V  ,  UQ2V  ,  UP2V  ,  UE2V  , 

&  UW  ,  UNV  ,  UTV  ,  UQV  ,  UPV  ,  UEV  ,  UAV  , 

&  UV2N  ,  UT2N  ,  UQ2N  ,  UP2N  ,  UE2N  , 

&  UNN  ,  UVN  ,  UTN  ,  UQN  ,  UPN  ,  UEN  ,  UAN  , 

&  UV2T  ,  UN2T  ,  UQ2T  ,  UP2T  ,  UE2T  , 

&  UTT  ,  UVT  ,  UNT  ,  UQT ,  UPT  ,  UET  ,  UAT  , 

&  UV2Q ,  UN2Q ,  UT2Q ,  UP2Q ,  UE2Q , 

&  UQQ ,  UVQ ,  UNQ ,  UTQ ,  UPQ ,  UEQ ,  UAQ , 

&  UV2P ,  UN2P ,  UT2P ,  UQ2P  ,  UE2P , 

&  UPP ,  UVP ,  UNP ,  UTP ,  UQP ,  UEP  ,  UAP 
c 

c  READ  INPUT  VARIABLES 
c 

open(ll/file='UA3.in',status='old') 

c 

read(ll,*)  TTxp 

read(ll/Jt')  date 

read  (11,*)  model 
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read(ll,*)  prop 

read(ll,*)  EHPxp 

read(ll,*)  SHPxp 

read(ll,*)  comment 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  ro  ,  nu 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  Uro,  Unu 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  L ,  S ,  D 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(U,*)  UL  ,  US  ,  UD 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  VCRknot ,  RTCR 

read(ll,*) 

read(ll,*) 

read(ll,*)  URTCR 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  Vknot ,  RT  ,  DelR 

read(ll,*) 

read(ll,*) 

read(ll,*)  nrpm  ,  T  ,  Qinlb 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  UV  ,  URT  ,  UDelR ,  Un  ,  UT  ,  UQ 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  JVCT  ,  tdCT ,  etaCT 
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read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  nusea 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*)  Uros  ,  Unus 
read(ll,*) 
read(ll,*) 
read(ll,*) 

read(ll,*)  LS  ,  VSknot 
read(ll,*) 
read(ll,*) 
read(ll,*) 

read  (11,*)  ULS  ,  USS  ,  UDS 

read(ll,*) 

read(ll,*) 

read(ll,*) 

read(ll,*j  UVfsm  ,  UNfsm  ,  UTfsm  ,  UQfsm  ,  UPfsm 

read(ll,*) 

read(ll,*) 

read  (11,*) 

read(ll,*)  CA ,  dCA 
c 

close(ll,status='keep') 

c 

c  About  notation  : 

c  U  stands  for  relative  Uncertainty 

c  fsm  stands  for  Full-Scale  Measurement  uncertainty 
c 

c  PRELIMINARY  TRANSFORMATIONS 

c 

c  Rescale  nu  and  nusea 

c 

nu  =  nu  /  100000. 
nusea  =  nusea  /  100000. 
c 

c  Transform  speeds  from  knots  to  ft/sec 

c 

VCR  =  1.6878  *  VCRknot 
V  =  1.6878  *  Vknot 
VS  =  1.6878  *  VSknot 
c 


36 


c  Transform  rpm  into  rps 

c 

n  =  nrpm  /  60. 
c 

c  Transform  torque  from  in-lb  to  ft-lb 
c 

Q  =  Qinlb  /  12. 
c 

c  Transform  input  percent  uncertainties 

c 

Uro  =  0.01  *  Uro 
Unu  =  0.01  *  Unu 
c 

UL  =  0.01  *  UL 
US  =  0.01  *  US 
UD  =  0.01  *  UD 
c 

URTCR  =  0.01  *  URTCR 
UV  =  0.01  *  uv 
URT  =  0.01  *  URT 
UDelR  =  0.01  *  UDelR 
c 

Un  =  0.01  *  Un 
UT  =  0.01  *  UT 
UQ  =  0.01  *  UQ 
c 

Uros  =  0.01  *  Uros 
Unus  =  0.01  *  Unus 
c 

ULS  =  0.01  *  ULS 
USS  =  0.01  *  USS 
UDS  =  0.01  *  UDS 
c 

UVfsm  =  0.01  *  UVfsm 
UNfsm  =  0.01  *  UNfsm 
UTfsm  =  0.01  *  UTfsm 
UQfsm  =  0.01  *  UQfsm 
UPfsm  =  0.01  *  UPfsm 
c 

c  MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS 
c 

c  Compute  CTCR  ,  CFCR  &  CR 
c 

CTCR  =  2.  *  RTCR  /  (  ro  *  S  *  VCR  *  VCR  ) 

CFCR  =  0.075  /  ( LOG10(  VCR  *  L  /  nu  )  -  2.  )**2 
CR  =  CTCR  -  CFCR 
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MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


c 


c 


c 


c 

c 

c 

c 

c 


c 

c 

c 


Compute  Ri ,  CT  &  CF 


Ri  =  RT  -  DelR 

CT  =  2.  *Ri/(ro*S*V*V) 

CF  =  0.075  /  (  LOG10(  V  *  L  /  nu  )  -  2.  )**2 

Compute  JV ,  td=l-t  &  eta=etaD 

JV  =  V/  (n*D) 
td  =  Ri  /  T 

eta  =  V  *  Ri  /  (  6.2831853  *n*Q) 

PRELIMINARY  CALCULATIONS  FOR  MODEL-SCALE  UNCERTAINTIES 


Uro2  =  Uro  *  Uro 
Unu2  =  Unu  *  Unu 

UL2  =  UL  *  UL 
US2  =  US  *  US 
UD2  =  UD  *  UD 

URTCR2  =  URTCR  *  URTCR 
UV2  =  UV  *  UV 

URi  =  (  URT  *  RT  )**2  +  (  UDelR  *  DelR  )**2 
URi  =  SQRT(  URi )  /  Ri 
URi2  =  URi  *  URi 


Un2  =  Un  *  Un 
UT2  =  UT  *  UT 
UQ2  =  UQ  *  UQ 

MODEL-SCALE  UNCERTAINTIES  :  RESISTANCE  TESTS 

Compute  dCTCR  /  CTCR 

UCTCR2  =  URTCR2  +  Uro2  +  US2  +  4.  *  UV2 
UCTCR  =  100.  *  SQRT(  UCTCR2  ) 

Compute  dCFCR  /  CFCR 

UCFCR2  =  10.  *  CFCR  *  (  UV2  +  UL2  +  Unu2  ) 

UCFCR  =  100.  *  SQRT(  UCFCR2  ) 
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c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


Compute  dCRA2  &  dCR  /  CR 

dCR2  =  UV2  *  (  2.  *  CTCR  -  CFCR  *  SQRT(  10.  *  CFCR  )  )**2 
dCR2  =  dCR2  +  (  URTCR2  +  Uro2  +  US2  )  *  CTCR**2 
dCR2  =  dCR2  +  (  UL2  +  Unu2  )  *  10.  *  CFCR**3 
UCR  =  100.  *  SQRT(  dCR2  )  /  CR 

MODEL-SCALE  UNCERTAINTIES  :  PROPULSION  TESTS 


Compute  dCT  /  CT 

UCT2  =  URi2  +  Uro2  +  US2  +  4.  *  UV2 
UCT  =  100.  *  SQRT(  UCT2  ) 


Compute  dCF  /  CF 

UCF2  =  10.  *  CF  *  (  UV2  +  UL2  +  Unu2  ) 
UCF  =  100.  *  SQRT(  UCF2  ) 


Compute  dJV  /  JV ,  dtd  /  td  &  deta  /  eta 

UJV2  =  UV2  +  Un2  +  UD2 

Utd2  =  URi2  +  UT2 

Ueta2  =  UV2  +  URi2  +  Un2  +  UQ2 

UJV  =  100.  *  SQRT(  UJV2  ) 

Utd  =  100.  *  SQRT(  Utd2  ) 

Ueta  =  100.  *  SQRT(  Ueta2  ) 

PRELIMINARY  CALCULATIONS  FOR  FULL-SCALE  UNCERTAINTIES 
Compute  CFS  &  CTS 

CFS  =  0.075  /  (  LOG10(  VS  *  LS  /  nusea  )  -  2.  )**2 
CIS  =  CFS  +  CR  +  CA 

Compute  (drhosea / rhosea) A2+ (dSS / SS) A2  &  (dJV/JV)A2+(dDS/DS)A2 


UroS2  =  Uros**2  +  USS**2 
UJD2  =  UJV2  +  UDS**2 

Compute  UCLnu2  &  Gamma 

cofCFS  =  10.  *  CFS**3 

UCLnu2  =  (  dCR2  +  dCA**2  )  /  cofCFS  +  ULS**2  +  Unus**2 
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cofCFS  =  SQRT(  cofCFS  ) 

Gamma  =  cofCFS  /  CTS 

Compute  sigmaj  ,  sigmat  &  sigmeta 

sigmaj  =  cofCFS  *  JVCT  /  JV 
sigmat  =  cofCFS  *  tdCT  /  td 
sigmeta  =  cofCFS  *  etaCT  /  eta 

Compute  squares  of  full-scale  measurement  uncertainties 

UVfsm2  =  UVfsm  *  UVfsm 
UNfsm2  =  UNfsm  *  UNfsm 
UTfsm2  =  UTfsm  *  UTfsm 
UQfsm2  =  UQfsm  *  UQfsm 
UPfsm2  =  UPfsm  *  UPfsm 

FULL-SCALE  UNCERTAINTIES  @  a  GIVEN  SPEED 

UN2V  =  UVfsm2  *  (  1.  +  sigmaj  )**2  +  UCLnu2  *  sigmaj**2 
UN2V  =  UN2V  +  UJD2  +  UNfsm2 

UT2V  =  UVfsm2  *  (  2.  -  Gamma  +  sigmat  )**2 
UT2V  =  UT2V  +  UCLnu2  *  (  Gamma  -  sigmat  )**2 
UT2V  =  UT2V  +  UroS2  +  Utd2  +  UTfsm2 


UQ2V  —  UVfsm2  *  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2 
UQ2V  =  UQ2V  +  UCLnu2  *  (  Gamma  -  sigmeta  +  sigmaj  )**2 
UQ2V  =  UQ2V  +  UroS2  +  Ueta2  +  UJD2  +  UQfsm2 

UP2V  =  UVfsm2  *  (  3.  -  Gamma  +  sigmeta  )**2 
UP2V  =  UP2V  +  UCLnu2  *  (  Gamma  -  sigmeta  )**2 
UP2V  =  UP2V  +  UroS2  +  Ueta2  +  UPfsm2 

UE2V  =  UVfsm2  *  (  3.  -  Gamma  )**2  +  UCLnu2  *  Gammaw2  +  UroS2 

UW  =  100.  *  UVfsm 
UNV  =  100.  *  SQRT(  UN2V  ) 

UTV  =  100.  *  SQRT(  UT2V  ) 

UQV  =  100.  *  SQRT(  UQ2V  ) 

UPV  =  100.  *  SQRT(  UP2V  ) 

UEV  =  100.  *  SQRT(  UE2V  ) 

UAV  =  20.  *  SQRT(  UVfsm2  +  UN2V  +  UT2V  +  UQ2V  +  UP2V  ) 
FULL-SCALE  UNCERTAINTIES  @  a  GIVEN  rpm 


UV2N  =  UNfsm2  +  UCLnu2  *  sigmaj**2  +  UJD2 
UV2N  =  UV2N  /  (  1.  +  sigmaj  )**2  +  UVfsm2 
c 

UT2N  =  (  UNfsm2  +  UJD2  )  *  (  2.  -  Gamma  +  sigmat  )**2 
UT2N  =  UT2N  +  UCLnu2  *  (  Gamma  -  sigmat  +  2.  *  sigmaj  )**2 
UT2N  =  UT2N  /  (  1.  +  sigmaj  )**2  +  UroS2  +  Utd2  +  UTfsm2 
c 

UQ2N  =  UNfsm2  *  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2 
UQ2N  =  UQ2N  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**2 
UQ2N  =  UQ2N  +  UJD2  *  (  3.  -  Gamma  +  sigmeta  )**2 
UQ2N  =  UQ2N  /  ( 1.  +  sigmaj  )**2  +  UroS2  +  Ueta2  +  UQfsm2 
c 

UP2N  =  (  UNfsm2  +  UJD2  )  *  (  3.  -  Gamma  +  sigmeta  )**2 
UP2N  =  UP2N  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**2 
UP2N  =  UP2N  /  ( 1.  +  sigmaj  )**2  +  UroS2  +  Ueta2  +  UPfsm2 
c 

UE2N  =  (  UNfsm2  +  UJD2  )  *  (  3.  -  Gamma  )**2 
UE2N  =  UE2N  +  UCLnu2  *  (  Gamma  +  3.  *  sigmaj  )**2 
UE2N  =  UE2N  /  ( 1.  +  sigmaj  )**2  +  UroS2 
c 

UNN  =  100.  *  UNfsm 
UVN  =  100.  *  SQRT(  UV2N  ) 

UTN  =  100.  *  SQRT(  UT2N  ) 

UQN  =  100.  *  SQRT(  UQ2N  ) 

UPN  =  100.  *  SQRT(  UP2N  ) 

UEN  =  100.  *  SQRT(  UE2N  ) 

UAN  =  20.  *  SQRT(  UV2N  +  UNfsm2  +  UT2N  +  UQ2N  +  UP2N  ) 
c 

c  FULL-SCALE  UNCERTAINTIES  @  a  GIVEN  THRUST 
c 

UV2T  =  UTfsm2  +  UCLnu2  *  (  Gamma  -  sigmat  )**2  +  UroS2  +  Utd2 
UV2T  =  UV2T  /  (  2.  -  Gamma  +  sigmat  )**2  +  UVfsm2 
c 

UN2T  =  (  UTfsm2  +  UroS2  +  Utd2  )  *  ( 1.  +  sigmaj  )**2 
UN2T  =  UN2T  +  UCLnu2  *  (  Gamma  -  sigmat  +  2.  *  sigmaj  )**2 
UN2T  =  UN2T  /  (  2.  -  Gamma  +  sigmat  )**2  +  UJD2  +  UNfsm2 
c 

UQ2T  =  (  UTfsm2  +  Utd2  )  *  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2 
UQ2T  =  UQ2T  +  (  4.  *  UCLnu2  +  UroS2  )  *  (  sigmat  -  sigmeta  +  sigmaj  )**2 
UQ2T  =  UQ2T  /  (  2.  -  Gamma  +  sigmat  )**2  +  UJD2  +  Ueta2  +  UQfsm2 
c 

UP2T  =  (  UTfsm2  +  Utd2  )  *  (  3.  -  Gamma  +  sigmeta  )**2 

UP2T  =  UP2T  +  UCLnu2  *  (  Gamma  +  2.  *  sigmeta  -  3.  *  sigmat  )**2 

UP2T  =  UP2T  +  UroS2  *  ( 1.  +  sigmeta  -  sigmat  )**2 

UP2T  =  UP2T  /  (  2.  -  Gamma  +  sigmat  )**2  +  Ueta2  +  UPfsm2 
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UE2T  =  (  UTfsm2  +  Utd2  )  *  (  3.  -  Gamma  )**2 
UE2T  =  UE2T  +  UCLnu2  *  (  Gamma  -  3.  *  sigmat  )**2 
UE2T  =  UE2T  +  UroS2  *  ( 1.  -  sigmat  )**2 
UE2T  =  UE2T  /  (  2.  -  Gamma  +  sigmat  )**2 

UTT  =  100.  *  UTfsm 
UVT  =  100.  *  SQRT(  UV2T  ) 

UNT  =  100.  *  SQRT(  UN2T  ) 

UQT  =  100.  *  SQRT(  UQ2T  ) 

UPT  =  100.  *  SQRT(  UP2T  ) 

UET  =  100.  *  SQRT(  UE2T  ) 

UAT  =  20.  *  SQRT(  UV2T  +  UN2T  +  UTfsm2  +  UQ2T  +  UP2T  ) 

FULL-SCALE  UNCERTAINTIES  @  a  GIVEN  TORQUE 

UV2Q  =  UQfsm2  +  UCLnu2  *  (  Gamma  -  sigmeta  +  sigmaj  )**2 
UV2Q  =  UV2Q  +  UroS2  +  Ueta2  +  UJD2 

UV2Q  =  UV2Q  /  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2  +  UVfsm2 

UN2Q  =  (  UQfsm2  +  UroS2  +  Ueta2  )  *  ( 1.  +  sigmaj  )**2 

UN2Q  =  UN2Q  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**2 

UN2Q  =  UN2Q  +  UJD2  *  (  3.  -  Gamma  +  sigmeta  )**2 

UN2Q  =  UN2Q  /  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2  +  UNfsm2 

UT2Q  =  (  UQfsm2  +  Ueta2  +  UJD2  )  *  (  2.  -  Gamma  +  sigmat  )**2 
UT2Q  =  UT2Q  +  (  4.  *  UCLnu2  +  UroS2  )  *  (  sigmat  -  sigmeta  +  sigmaj  )**2 
UT2Q  =  UT2Q  /  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2  +  Utd2  +  UTfsm2 

UP2Q  =  •(  UQfsm2  +  UJD2  )  *  (  3.  -  Gamma  +  sigmeta  )**2 
UP2Q  =  UP2Q  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**2 
UP2Q  =  UP2Q  +  (  UroS2  +  Ueta2  )  *  ( 1.  +  sigmaj  )**2 
UP2Q  =  UP2Q  /  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2  +  UPfsm2 

UE2Q  =  (  UQfsm2  +  UJD2  )  *  (  3.  -  Gamma  )**2 
UE2Q  =  UE2Q  +  UCLnu2  *  (  Gamma  +  3.  *  sigmaj  )**2 
UE2Q  =  UE2Q  +  UroS2  *  ( 1.  +  sigmaj  )**2 
UE2Q  =  UE2Q  /  ( 2.  -  Gamma  -  sigmaj  )**2 

UQQ  =  100.  *  UQfsm 
UVQ  =  100.  *  SQRT(  UV2Q  ) 

UNQ  =  100.  *  SQRT(  UN2Q  ) 

UTQ  =  100.  *  SQRT(  UT2Q  ) 

UPQ  =  100.  *  SQRT(  UP2Q ) 

UEQ  =  100.  *  SQRT(  UE2Q ) 

UAQ  =  20.  *  SQRT(  UV2Q  +  UN2Q  +  UT2Q  +  UQfsm2  +  UP2Q  ) 


FULL-SCALE  UNCERTAINTIES  @  a  GIVEN  SHAFT  HORSEPOWER 


UV2P  =  UPfsm2  +  UCLnu2  *  (Gamma-sigmeta)**2  +  UroS2  +  Ueta2 
UV2P  =  UV2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  UVfsm2 

UN2P  =  (  UPfsm2  +  UroS2  +  Ueta2  )  *  ( 1.  +  sigmaj  )**2 
UN2P  =  UN2P  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**2 
UN2P  =  UN2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  UJD2  +  UNfsm2 

UT2P  =  (  UPfsm2  +  Ueta2  )  *  (  2.  -  Gamma  +  sigmat  )**2 
UT2P  =  UT2P  +  UCLnu2  *  (Gamma+2.*sigmeta-3.*sigmat)**2 
UT2P  =  UT2P  +  UroS2  *  ( 1.  +  sigmeta  -  sigmat  )**2 
UT2P  =  UT2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  Utd2  +  UTfsm2 

UQ2P  =  UPfsm2  *  (  2.  -  Gamma  +  sigmeta  -  sigmaj  )**2 
UQ2P  =  UQ2P  +  UCLnu2  *  (  Gamma  -  sigmeta  +  3.  *  sigmaj  )**2 
UQ2P  =  UQ2P  +  (  UroS2  +  Ueta2  )  *  ( 1.  +  sigmaj  )**2 
UQ2P  =  UQ2P  /  (  3.  -  Gamma  +  sigmeta  )**2  +  UJD2  +  UQfsm2 

UE2P  =  (  UPfsm2  +  Ueta2  )  *  (  3.  -  Gamma  )**2 
UE2P  =  UE2P  +  (  9.  *  UCLnu2  +  UroS2  )  *  sigmeta**2 
UE2P  =  UE2P  /  (  3.  -  Gamma  +  sigmeta  )**2 

UPP  =  100.  *  UPfsm 
UVP  =  100.  *  SQRT(  UV2P  ) 

UNP  =  100.  *  SQRT(  UN2P  ) 

UTP  =  100.  *  SQRT(  UT2P  ) 

UQP  =  100.  *  SQRT(  UQ2P  ) 

UEP  =  100.  *  SQRT(  UE2P ) 

UAP  =  20.  *  SQRT(  UV2P  +  UN2P  +  UT2P  +  UQ2P  +  UPfsm2  ) 

WRITE  INPUT  VARIABLES  &  OUTPUT  RESULTS 

Express  relative  uncertainties  in  percent 

Uro  =  100.  *  Uro 
Unu  =  100.  *  Unu 

UL  =  100.  *  UL 
US  =  100.  *  US 
UD  =  100.  *  UD 

URTCR  =  100.  *  URTCR 
UV  =  100.  *  uv 
URT  =  100.  *  URT 
UDelR  =  100.  *  UDelR 


Un  =  100.  *  Un 
UT  =  100.  *  UT 
UQ  =  100.  *  UQ 

Uros  =  100.  *  Uros 
Unus  =  100.  *  Unus 

ULS  =  100.  *  ULS 
USS  =  100.  *  USS 
UDS  =  100.  *  UDS 

UVfsm  =  100.  *  UVfsm 
UNfsrn  =  100.  *  UNfsm 
UTfsm  =  100.  *  UTfsm 
UQfsm  =  100.  *  UQfsm 
UPfsm  =  100.  *  UPfsm 

UCA  =  100.  *  dCA  /  CA 

open(12,file='UA3.out',status='new') 

write(12,*)  TTxp 
write(12,*)  date 
write(12,*)  model 
write(12,*)  prop 
write(12,*)  EHPxp 
write(12,*)  SHPxp 
write(12,*)  comment 

write(12,*) 

write(12,*)  ’  INPUT  VARIABLES' 
write(12,*) 

write(12,*) 

write(12,*)  '  TANK-WATER  PROPERTIES 

write(12,*) 

write(12,101)  ro 

write(12,102)  nu 

write(12,103)  Uro 

write(12,104)  Unu 

write(12,*) 

write(12,*)  '  MODEL  GEOMETRY' 

write(12,*) 

write(12,105)  L 
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write(12/106)  S 
write(12,107)  D 
write(12/108)  UL 
write(12/109)  US 
write(12,110)  UD 

write(12,*) 

write(12,*)  '  MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS' 
write(12,*) 

write(12,lll)  VCRknot 
write(12,112)  RTCR 
write(12,U3)  URTCR 

write(12,*) 

write(12,*)  '  MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS' 

write  (12,*) 

write(12,114)  Vknot 

write(12,115)  RT 

write(12,116)  DelR 

write(12,117)  nrpm 

write(12,118)  T 

write(12,119)  Qinlb 

write(12,*) 

write(12,*)  ’  UNCERTAINTIES  OF  MODEL-SCALE  MEASUREMENTS’ 

write(12,*) 

write(12/120)  UV 

write(12,121)  URT 

write(12/122)  UDelR 

write(12,123)  Un 

write(12/124)  UT 

write(12/125)  UQ 

write(12,*) 

write(12/*)  ’  OTHER  MODEL-SCALE  VARIABLES' 

write(12,*) 

write(12/126)  JVCT 

write(12,127)  tdCT 

write(12,128)  etaCT 

write(12,*) 

write(12,*)  '  SEA-WATER  PROPERTIES' 
write  (12,*) 
write(12,151)  nusea 
write(12/152)  Uros 
write(12,153)  Unus 
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write(12,*) 

write(12,*)  ’  FULL-SCALE  SHIP' 

write(12,*) 

write(12,154)  LS 

write(12/155)  VSknot 

write(12,156)  ULS 

write(12/157)  USS 

write  (12, 158)  UDS 

write(12,*) 

write(12,*)  ’  UNCERTAINTIES  OF  FULL-SCALE  MEASUREMENTS 

write(12,*) 

write(12,159)  UVfsm 

write(12,160)  UNfsm 

write(12,161)  UTfsm 

write(12,162)  UQfsm 

write(12,163)  UPfsm 

write(12,*) 

write(12,*)  ’  SCALING  ALLOWANCE' 
write(12,*) 
write(12,171)  CA 
write(12,172)  dCA 

write(12,*) 

write(12,*) 

write(12,*)  '  OUTPUT  VARIABLES' 
write(12,*) 

write(12,*) 

write(12,*)  ’  MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS' 

write(12,*) 

write(12,201)  CTCR 

write(12,202)  CFCR 

write( 12,203)  CR 

write(12,*) 

write(12,)f)  ’  MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS' 

write(12,*) 

write(12,204)  CT 

write(12,205)  CF 

write(12,206)  CR 

write(12,If) 

write(12,207)  JV 

write(12,208)  td 
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write(12,209)  eta 
write(12,*) 

write (12,*)  '  UNCERTAINTIES  OF  MODEL-SCALE  VARIABLES  : 

RESISTANCE  TESTS' 

write(12,*) 

write(12,210)  UCTCR 
write(12,211)  UCFCR 
write(12,212)  UCR 

write(12,*) 

write(12,*)  '  UNCERTAINTIES  OF  MODEL-SCALE  VARIABLES  : 

PROPULSION  TESTS' 

write(12,*) 

write(12,213)  UCT 

write(12,214)  UCF 

write(12,215)  UCR 

write(12,216)  UCA 

write(12,*) 

write(12,217)  UJV 

write(12,218)  Utd 

write(12,219)  Ueta 

write(12,*) 

write(12,*)  '  FULL-SCALE  VARIABLES' 

write(12,*) 

write( 12,220)  CTS 

write(12,221)  CFS 

write(12,222)  CR 

write(12,223)  CA 

write(12,*) 

write(12,*)  '  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 

SPEED' 

write(12,*) 

write(12,311)  UW 

write(12,312)  UNV 

write(12,313)  UTV 

write(12,314)  UQV 

write(12,315)  UPV 

write(12,316)  UEV 

write(12,317)  UAV 

write(12,*) 

write(12,*)  '  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 
rpm' 
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write(12,*) 
write(12,321)  UVN 
write(12,322)  UNN 
write(12,323)  UTN 
write(12,324)  UQN 
write(12,325)  UPN 
write(12,326)  UEN 
write(12,327)  UAN 

write(12,*) 

write(12,*)  ’  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 

THRUST’ 

write(12,*) 

write(12,331)  UVT 

write(12,332)  UNT 

write( 12,333)  UTT 

write(12,334)  UQT 

write(12,335)  UPT 

write(12,336)  UET 

write(12,337)  UAT 

write(12,*) 

write(12,*)  ’  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 

TORQUE' 

write(12,*) 

write(12,341)  UVQ 

write(12,342)  UNQ 

write(12,343)  UTQ 

write(12,344)  UQQ 

write(12,345)  UPQ 

write( 12,346)  UEQ 

write(12,347)  UAQ 

write  (12,*) 

write(12,*)  ’  UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN 
SHP' 

write  (12,*) 
write(12,351)  UVP 
write(12,352)  UNP 
write(12,353)  UTP 
write( 12,354)  UQP 
write(12,355)  UPP 
write(12,356)  UEP 
write(12,357)  UAP 

close(12,status=’keep') 
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FORMATS 


c 
c 
c 

101  format('  water  density  (slug/ft**3) :  ',F11.3) 

102  format('  kinematic  viscosity  coefficient  (ft**2/sec)  :  ’,E11.4) 

103  format('  percent  uncertainty  of  density  :  ',F8.2) 

104  format('  percent  uncertainty  of  viscosity  :  ',F8.2) 
c 

105  format('  length  of  model  (ft) :  ’,F8.3) 

106  format('  wetted  area  (ft**2) :  ',F8.3) 

107  format('  diameter  of  propeller  (ft) :  ’,F8.4) 

108  format('  percent  uncertainty  of  length  :  ',F8.2) 

109  format('  percent  uncertainty  of  wetted  surface  :  ’,F8.2) 

110  format(’  percent  uncertainty  of  prop  diameter  :  ',F8.2) 
c 

111  format('  carriage  speed  (knots) :  ’,F8.2) 

112  format(’  drag  (lbs) :  ’,F8.2) 

113  format('  percent  uncertainty  of  drag  :  ',F8.2) 

c 

114  format('  carriage  speed  (knots) :  ’,F8.2) 

115  format('  drag  RT  (lbs) :  ’,F8.2) 

116  formate  drag  DeltaR  (lbs) :  ’JF8.2) 

117  format('  propeller  rpm  :  ',F8.2) 

118  format('  propeller  thrust  (lbs) :  ’,F8.2) 

119  format('  propeller  torque  (in-lbs)  :  ',F8.2) 
c 

120  format('  percent  uncertainty  of  model  speed  :  ',F8.2) 

121  format('  percent  uncertainty  of  drag  RT  :  ',F8.2) 

122  formate  percent  uncertainty  of  drag  DeltaR  :  ’,F8.2) 

123  format('  percent  uncertainty  of  propeller  rpm  :  ’,F8.2) 

124  format('  percent  uncertainty  of  prop  thrust :  ',F8.2) 

125  format('  percent  uncertainty  of  prop  torque  :  ',F8.2) 

c 

126  formate  slope  d  JV  /  d  CT  :  ',F9.3) 

127  format(’  slope  d  (1-t)  /  d  CT  :  ',F9.3) 

128  format('  slope  d  etaD  /  d  CT  :  ’,F9.3) 
c 

151  format('  kinematic  viscosity  coefficient  (ft**2/sec)  :  ’,E11.4) 

152  format('  percent  uncertainty  of  density  :  ’,F8.2) 

153  format('  percent  uncertainty  of  viscosity  :  ’,F8.2) 
c 

154  format('  ship  length  (ft) :  ’,F8.2) 

155  format('  ship  speed  (knots) :  '.Fill) 

c 

156  format('  percent  uncertainty  of  ship  length  :  ',F8.2) 

157  format('  percent  uncertainty  of  wetted  surface  :  ',F8.2) 
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158  format('  percent  uncertainty  of  prop  diameter  :  ',F8.2) 
c 

159  format('  percent  uncertainty  of  ship  speed  :  ',F8.2) 

160  formate  percent  uncertainty  of  prop  rpm  :  ',F8.2) 

161  format('  percent  uncertainty  of  prop  thrust :  ',F8.2) 

162  format('  percent  uncertainty  of  prop  torque  :  ’,F8.2) 

163  format('  percent  uncertainty  of  SHP  :  ',F8.2) 

c 

171  format('  correlation  allowance  :  ',F11.5) 

172  format('  uncertainty  of  allowance  :  ’,F11.5) 
c 

201  format('  total  resistance  coefficient  CT  :  ’,F10.5) 

202  format(’  friction  resistance  coefficient  CF  :  ',F10.5) 

203  format('  residuary  resistance  coefficient  CR  :  ',F10.5) 

c 

204  format('  total  resistance  coefficient  CT  :  ’,F10.5) 

205  formate  friction  resistance  coefficient  CF  :  ',F10.5) 

206  format(’  residuary  resistance  coefficient  CR  :  ’,F10.5) 

c 

207  format('  advance  ratio  JV  :  ',F8.2) 

208  format('  thrust-deduction  factor  1-t :  ’,F8.2) 

209  format('  propulsive  efficiency  etaD  :  ’,F8.2) 
c 

210  format('  percent  uncertainty  of  CT  :  ',F8.2) 

211  formate  percent  uncertainty  of  CF  :  ’,F8.2) 

212  formate  percent  uncertainty  of  CR  :  ’,F8.2) 
c 

213  format('  percent  uncertainty  of  CT  :  ',F8.2) 

214  formate  percent  uncertainty  of  CF  :  ’,F8.2) 

215  formate  percent  uncertainty  of  CR  :  ',F8.2) 

216  format('  percent  uncertainty  of  CA  :  ’,F8.2) 
c 

217  format('  percent  uncertainty  of  JV  :  ’,F8.2) 

218  formate  percent  uncertainty  of  1-t :  ',F8.2) 

219  formate  percent  uncertainty  of  etaD  :  ’,F8.2) 
c 

220  format('  total  resistance  coefficient  CT  :  ',F10.5) 

221  format('  friction  resistance  coefficient  CF  :  ’,F10.5) 

222  format('  residuary  resistance  coefficient  CR  :  ’,F10.5) 

223  format('  correlation  allowance  coefficient  CA  :  ’,F10.5) 
c 

311  format('  percent  uncertainty  of  ship  speed  :  ’,F8.2) 

312  format('  percent  uncertainty  of  prop  rpm  :  ’,F8.2) 

313  format('  percent  uncertainty  of  prop  thrust :  ',F8.2) 

314  format('  percent  uncertainty  of  prop  torque  :  ’,F8.2) 

315  format('  percent  uncertainty  of  SHP  :  ’,F8.2) 
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316  format('  percent  uncertainty  of  EHP  :  ',F8.2) 

317  format('  percent  overall  uncertainty  :  ',F8.2) 

c 

321  format('  percent  uncertainty  of  ship  speed  :  ’,F8.2) 

322  format('  percent  uncertainty  of  prop  rpm  :  ’,F8.2) 

323  format('  percent  uncertainty  of  prop  thrust :  ’,F8.2) 

324  format('  percent  uncertainty  of  prop  torque  :  ’,F8.2) 

325  formate  percent  uncertainty  of  SHP  :  ’,F8.2) 

326  formate  percent  uncertainty  of  EF1P  :  ’,F8.2) 

327  format('  percent  overall  uncertainty  :  ',F8.2) 

c 

331  formate  percent  uncertainty  of  ship  speed  :  ',F8.2) 

332  format('  percent  uncertainty  of  prop  rpm  :  ’,F8.2) 

333  format('  percent  uncertainty  of  prop  thrust :  ',F8.2) 

334  format('  percent  uncertainty  of  prop  torque  :  ',F8.2) 

335  format('  percent  uncertainty  of  SHP  :  ',F8.2) 

336  formate  percent  uncertainty  of  EHP  :  ',F8.2) 

337  format('  percent  overall  uncertainty  :  ',F8.2) 

c 

341  format('  percent  uncertainty  of  ship  speed  :  ',F8.2) 

342  format('  percent  uncertainty  of  prop  rpm  :  ’,F8.2) 

343  format('  percent  uncertainty  of  prop  thrust :  ',F8.2) 

344  format('  percent  uncertainty  of  prop  torque  :  ',F8.2) 

345  format('  percent  uncertainty  of  SHP  :  ',F8.2) 

346  formate  percent  uncertainty  of  EHP  :  ',F8.2) 

347  format('  percent  overall  uncertainty  :  ',F8.2) 

c 

351  format('  percent  uncertainty  of  ship  speed  :  ’,F8.2) 

352  format('  percent  uncertainty  of  prop  rpm  :  ',F8.2) 

353  format('  percent  uncertainty  of  prop  thrust :  ’,F8.2) 

354  format('  percent  uncertainty  of  prop  torque  :  ',F8.2) 

355  format('  percent  uncertainty  of  SHP  :  ',F8.2) 

356  formate  percent  uncertainty  of  EHP  :  ',F8.2) 

357  format('  percent  overall  uncertainty  :  ',F8.2) 

c 

stop 

c 

end 
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EXAMPLE  INPUT  FILE 


'  TOW-TANK  EXP  ' 

'  DATE:  XXXX  ' 

’  MODEL  No.  XXXX  ' 

'  PROPELLER  No.  XXXX  ' 

'  EHP  EXPERIMENT  No.  XXXX  ' 

'  SHP  EXPERIMENT  No.  XXXX  ' 

’  COMMENTS:  CASE  MFC  ' 

MODEL-SCALE  VARIABLES  AND  UNCERTAINTIES 

Tank- water  density  (rho)  and  kinematic  viscosity  (nu) 
rho  (slug/ ft**3)  nu  X  10**5  (ft**2/sec) 

1.9367 ,  1.084 

Percent  relative  uncertainties  of  rho  &  nu 
density  kinematic  viscosity 
0.07,  2.1 

Model  geometry 

length  (ft)  area  (ft**2)  prop  diameter  (ft) 

22.697,  138.179,  0.9986 

Percent  relative  uncertainties  of  model  geometry 
length  area  prop  diameter 

0.14 ,  0.71 ,  0.071 

RESISTANCE  (EHP)  TESTS 

speed  (knots)  drag  (lbs) 

6.0 ,  46.66 

Percent  relative  uncertainty  of  drag 


PROPULSION  TESTS 

speed  (knots)  RT  (lbs)  DeltaR 

18.0  ,  392.64  ,  70.1 

rpm  thrust  (lbs)  torque  (in-lbs) 

923.4,  551.0,  1501.0 

Percent  relative  uncertainties 

speed  RT  DeltaR  rpm  thrust  torque 

0.21,  1.7,  1.7,  0.42,  0.71,  0.71 
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Other  model  variables  :  slopes  of  1-t ,  JV  &  etaD  versus  CT 
dJV/dCT  d(l-t)/dCT  detaD/dCT 

-0.249,  0.067,  -0.015 

FULL-SCALE  VARIABLES  AND  UNCERTAINTIES 

Sea-water  kinematic  viscosity 
nu  X  10**5  (ft**2/sec) 

1.282 

Percent  relative  uncertainties  of  sea-water  properties 
density  kinematic  viscosity 

1.0 ,  2.0 

Full-scale  variables 

ship  length  (ft)  ship  speed  (knots) 

380.0 ,  25.0 

Percent  relative  uncertainties  of  full-scale  geometry 
length  area  prop  diameter 

0.5 ,  1.0 ,  0.1 

Percent  relative  uncertainties  of  full-scale  measurements 
speed  rpm  thrust  torque  SHP 

0.6 ,  0.4 ,  3.0 ,  0.9 ,  0.9 

SCALING  ALLOWANCE 
CA  dCA 

0.0005 ,  0.0001 
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EXAMPLE  OUTPUT  FILE 


TOW-TANK  EXP 
DATE:  XXXX 
MODEL  No.  XXXX 
PROPELLER  No.  XXXX 
EHP  EXPERIMENT  No.  XXXX 
SHP  EXPERIMENT  No.  XXXX 
COMMENTS:  CASE  MFC 


INPUT  VARIABLES 


TANK-WATER  PROPERTIES 

water  density  (slug/ft**3) :  1.937 

kinematic  viscosity  coefficient  (ft**2/sec) :  0.1084E-04 
percent  uncertainty  of  density  :  0.07 

percent  uncertainty  of  viscosity  :  2.10 

MODEL  GEOMETRY 

length  of  model  (ft) :  22.697 

wetted  area  (ft**2) :  138.179 

diameter  of  propeller  (ft) :  0.9986 

percent  uncertainty  of  length  :  0.14 

percent  uncertainty  of  wetted  surface  :  0.71 

percent  uncertainty  of  prop  diameter  :  0.07 

MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS 

carriage  speed  (knots) :  6.00 

drag  (lbs) :  46.66 

percent  uncertainty  of  drag :  1.70 

MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS 

carriage  speed  (knots) :  18.00 

drag  RT  (lbs) :  392.64 

drag  DeltaR  (lbs) :  70.10 

propeller  rpm :  923.40 

propeller  thrust  (lbs) :  551.00 

propeller  torque  (in-lbs) :  1501.00 

UNCERTAINTIES  OF  MODEL-SCALE  MEASUREMENTS 

percent  uncertainty  of  model  speed  :  0.21 
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percent  uncertainty  of  drag  RT :  1.70 

percent  uncertainty  of  drag  DeltaR :  1.70 

percent  uncertainty  of  propeller  rpm  :  0.42 

percent  uncertainty  of  prop  thrust :  0.71 

percent  uncertainty  of  prop  torque  :  0.71 

OTHER  MODEL-SCALE  VARIABLES 

slope  d  JV  /  d  CT :  -0.249 

slope  d  (1-t)  /  d  CT :  0.067 

slope  d  etaD  /  d  CT :  -0.015 

SEA- WATER  PROPERTIES 

kinematic  viscosity  coefficient  (ft**2/ sec)  :  0.1282E-04 
percent  uncertainty  of  density  :  1.00 

percent  uncertainty  of  viscosity  :  2.00 

FULL-SCALE  SHIP 

ship  length  (ft) :  380.00 

ship  speed  (knots) :  25.0 

percent  uncertainty  of  ship  length  :  0.50 

percent  uncertainty  of  wetted  surface  :  1.00 

percent  uncertainty  of  prop  diameter  :  0.10 

UNCERTAINTIES  OF  FULL-SCALE  MEASUREMENTS 

percent  uncertainty  of  ship  speed  :  0.60 

percent  uncertainty  of  prop  rpm  :  0.40 

percent  uncertainty  of  prop  thrust :  3.00 

percent  uncertainty  of  prop  torque  :  0.90 

percent  uncertainty  of  SHP  :  0.90 

SCALING  ALLOWANCE 

correlation  allowance  :  0.00050 

uncertainty  of  allowance  :  0.00010 


OUTPUT  VARIABLES 

MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS 

total  resistance  coefficient  CT  :  0.00340 

friction  resistance  coefficient  CF  :  0.00264 
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residuary  resistance  coefficient  CR :  0.00075 

MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS 

total  resistance  coefficient  CT  :  0.00261 

friction  resistance  coefficient  CF  :  0.00223 

residuary  resistance  coefficient  CR  :  0.00075 

advance  ratio  JV  :  1.98 

thrust-deduction  factor  1-t :  0.59 

propulsive  efficiency  etaD  :  0.81 

UNCERTAINTIES  OF  MODEL-SCALE  VARIABLES  :  RESISTANCE  TESTS 

percent  uncertainty  of  CT  :  1.89 

percent  uncertainty  of  CF  :  0.34 

percent  uncertainty  of  CR  :  8.58 

UNCERTAINTIES  OF  MODEL-SCALE  VARIABLES  :  PROPULSION  TESTS 

percent  uncertainty  of  CT  :  2.26 

percent  uncertainty  of  CF  :  0.32 

percent  uncertainty  of  CR  :  8.58 

percent  uncertainty  of  CA  :  20.00 

percent  uncertainty  of  JV  :  0.47 

percent  uncertainty  of  1-t :  2.22 

percent  uncertainty  of  etaD  :  2.27 

FULL-SCALE  VARIABLES 


total  resistance  coefficient  CT  :  0.00274 

friction  resistance  coefficient  CF  :  0.00149 

residuary  resistance  coefficient  CR  :  0.00075 

correlation  allowance  coefficient  CA  :  0.00050 

UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  SPEED 


percent  uncertainty  of  ship  speed  :  0.60 

percent  uncertainty  of  prop  rpm  :  0.87 

percent  uncertainty  of  prop  thrust :  6.01 

percent  uncertainty  of  prop  torque  :  5.33 

percent  uncertainty  of  SHP  :  5.47 

percent  uncertainty  of  EHP  :  4.90 

percent  overall  uncertainty  :  1.96 
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UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  rpm 


percent  uncertainty  of  ship  speed  :  0.87 

percent  uncertainty  of  prop  rpm  :  0.40 

percent  uncertainty  of  prop  thrust :  6.02 

percent  uncertainty  of  prop  torque  :  5.42 

percent  uncertainty  of  SHP  :  5.50 

percent  uncertainty  of  EHP  :  4.92 

percent  overall  uncertainty  :  1.97 

UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  THRUST 

percent  uncertainty  of  ship  speed  :  3.11 

percent  uncertainty  of  prop  rpm  :  3.11 

percent  uncertainty  of  prop  thrust :  3.00 

percent  uncertainty  of  prop  torque  :  4.48 

percent  uncertainty  of  SHP  :  6.60 

percent  uncertainty  of  EHP :  6.13 

percent  overall  uncertainty  :  1.92 

UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  TORQUE 

percent  uncertainty  of  ship  speed  :  2.76 

percent  uncertainty  of  prop  rpm  :  2.80 

percent  uncertainty  of  prop  thrust :  4.48 

percent  uncertainty  of  prop  torque  :  0.90 

percent  uncertainty  of  SHP :  3.19 

percent  uncertainty  of  EHP  :  2.82 

percent  overall  uncertainty  :  1.36 

UNCERTAINTIES  of  FULL-SCALE  PREDICTIONS  at  a  GIVEN  SHP 

percent  uncertainty  of  ship  speed  :  1.86 

percent  uncertainty  of  prop  rpm  :  1.87 

percent  uncertainty  of  prop  thrust :  4.35 

percent  uncertainty  of  prop  torque  :  2.10 

percent  uncertainty  of  SHP  :  0.90 

percent  uncertainty  of  EHP :  2.44 

percent  overall  uncertainty  :  1.12 
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APPENDIX  D: 

REPEATABILITY  OF  MODEL-SCALE  MEASUREMENTS 


The  precision  uncertainties  associated  with  submarine  model  tow-tank 
testing  are  investigated  by  considering  two  complete  tow-tank  test  series.  These  two 
test  series  are  representative  of  the  submarine  model  resistance  and  powering 
experimental  evaluations  currently  being  performed  in  the  NSWCCD  towing  tank. 
The  first  test  series  was  performed  from  01/2/96  to  2/20/96.  The  second  test  series 
was  performed  from  5/19/97  to  06/20/9 7.  The  current  model  test  speed  range  of  6-18 
knots  is  represented  in  these  tests,  as  well  as  the  current  data  collection 
instrumentation  and  calibration  techniques. 

The  precision  uncertainties  of  the  model  measurements  of  speed,  drag,  rpm, 
thrust,  and  torque  are  evaluated  two  ways.  First,  the  gage  calibrations  and 
instrumentation  were  analyzed  for  uncertainties.  The  drag,  thrust,  and  torque  gages 
are  calibrated  on  site  and  the  electronic  instrumentation  manufacturer  specifications 
are  examined.  These  precision  uncertainties  are  generally  very  small.  A  better 
assessment  of  the  model  measurement  precision  uncertainties  is  obtained  by 
analyzing  the  collected  test  data.  This  second  way  of  evaluating  the  measurement 
precision  uncertainties  takes  into  account  the  whole  data  collection  system, 
including  the  effects  of  changing  model  conditions  during  a  test  series,  water  flow 
variations,  variation  in  force  gages,  instrumentation  accuracy,  carriage  vibrations, 
and  computer  collection  and  recording  of  the  collected  model  speed,  model  drag, 
shaft  RPM,  shaft  thrust  and  shaft  torque  values.  The  methods  used  to  determine 
the  precision  uncertainties  of  the  model  provide  conservative  uncertainty  values 
for  use  in  the  global  uncertainty  analysis.  The  calibration  and  measurement 
precision  uncertainties  for  the  five  main  measured  model  quantities,  i.e.,  model 
speed  and  drag,  shaft  RPM,  thrust,  and  torque,  are  examined  below. 

Calibration  of  the  drag,  thrust,  and  torque  gages  is  completed  before  each  test 
series.  These  calibrations  consist  of  determining  a  unit/volt  factor  which  is  used  for 
converting  from  measured  volts  to  physical  units.  The  weights  used  for  the 
calibration  are  calibrated  every  two  years  to  a  tolerance  of  0.01%  of  the  nominal 
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weight  by  the  State  of  Maryland  Department  of  Agriculture  Weights  and  Measures 
Section.  The  tolerances  are  given  in  Table  2  through  5  of  the  NIST  Handbook  105-1 
(Revised  1990),  "Class  F  Tolerances  for  Field  Standards  Weights."  The  uncertainty 
of  a  single  measurement  within  a  calibration  is  calculated  by  first  taking  two  times 
the  standard  deviation  of  the  data  spots  within  that  calibration.  The  average  of 
these  standard  deviations  for  all  the  calibrations  is  calculated  next  and  represents  the 
precision  uncertainty  of  a  single  measurement  in  the  test  series.  This  also  provides 
verification  of  the  linearity  of  the  calibration  factor.  The  cal  factor  used  for  the  test  is 
determined  by  averaging  the  multiple  number  of  calibrations  for  each  gage.  The 
precision  error  associated  with  determining  the  cal  factor  is  given  by  two  times  the 
standard  deviation  of  the  set  of  averages  of  the  calibrations.  Figure  D.l.  presents 
calibration  examples  for  drag,  thrust,  and  torque.  The  instrumentation  currently 
being  used  for  the  model  test  measurements  is  presented  in  Table  D.l.  The 
calibration  precision  uncertainty  results  are  presented  in  Table  D.2.  and  D.3.  for  two 
test  series. 

A  typical  resistance  or  powering  experiment  will  consist  of  a  minimum  of 
approximately  20  data  spots  with  each  data  spot  consisting  of  the  average  value  of  5 
seconds  of  data  collection  at  400  samples/second  for  a  total  sample  of  2000  for  a  set 
speed,  drag,  and  RPM  (for  powering)  condition. 

The  precision  uncertainty  for  model  speed  is  calculated  using  all  the  speed 
data  from  the  tests  within  the  two  series.  Two  times  the  standard  deviation  of 
approximately  20  data  spots  for  each  test  was  calculated  to  determine  the  precision 
uncertainty  for  that  particular  test.  The  results  are  presented  in  Tables  D.4.  and  D.5. 
Model  speed  collected  for  a  typical  experiment  is  presented  in  Figure  D.2. 

The  precision  uncertainty  in  model  drag  is  calculated  by  using  measured 
values  of  drag  at  the  same  nominal  model  speed.  A  correction  is  applied  to  the 
measured  drag  to  reduce  the  effects  of  the  speed  variation  on  the  measured  drag 
uncertainty.  The  measured  drag  is  multiplied  by  (VN)2  /(VM)2,  where  VM  =  measured 
speed  and  VN  =  nominal  speed.  EHP  tests,  before  no-loads,  and  after  no-loads  are 
used  as  the  source  of  the  data.  This  drag  data  was  taken  throughout  the  test  series 
and  therefore  represents  a  large  data  set  for  the  uncertainty  analysis.  The  no-loads 
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are  collected  before  and  after  most  SHP  tests  to  find  the  model  drag  and  other  forces. 
No-loads  are  also  a  check  to  ensure  that  test  conditions  have  not  changed.  The 
precision  uncertainties  are  presented  in  Tables  D.6.  and  D.7.  Model  drag 
measurements  for  a  complete  test  series  are  presented  in  Figure  D.3. 

The  precision  uncertainties  for  model  thrust,  torque,  and  RPM  are 
determined  differently  from  the  model  drag  and  speed  uncertainties.  Typical 
submarine  powering  experiments  consist  of  varying  the  propeller  rpm  to  produce 
different  submarine  loadings.  The  precision  uncertainty  for  the  shaft  thrust,  torque, 
and  RPM  has  been  estimated  by  determining  the  variation  of  the  thrust,  torque,  and 
RPM  data  from  a  least  squares  curve  fit  through  the  data  spots  of  a  test.  Each  test 
contains  about  20  data  spots  which  comprise  a  range  of  thrust,  torque  and  RPM 
versus  total  drag  coefficient  (CT)  values.  The  thrust,  torque,  and  RPM  are  plotted 
against  CT  and  a  second-order  least  squared  curve  is  fit  through  each  set  of  data. 

Twice  the  standard  error  estimate  of  the  data  set  is  calculated  and  taken  to  be  the 
precision  uncertainty  for  that  test.  Examples  are  shown  in  Figures  D.4.,  D.5.,  and  D.6. 
The  results  are  presented  in  Tables  D.8.  and  D.9. 


The  following  table  summarizes  the  precision  uncertainties  for  tow-tank 
model  testing. 


Up  speed 

Up  drag 

Up  thrust 

Up  torque 

Up  rpm 

0.15% 

1.2% 

0.5% 

0.5% 

0.3% 

60 


,  Ibs/volt 


O  Drag  Ahead  Cal  Factor,  Ibs/volt 
■  —  Average 


Average  +2*Stdev 
Average  -2*Stdev 


Model  Speed,  ft/s 


Spot  # 

Fig.  D.2.  Example  of  model  speed  measurements  and  uncertainty 
for  one  experiment 
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Resi 


425  430  435  440  445  450  455  460  465 

Exp.  # 


Fig.  D.3.  Example  of  model  drag  measurements  and  uncertainty  for  a 
complete  test  series 
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Table  D.l.  Instrumentation  used  for  the  two  resistance  and  powering  tests 


INSTRUMENT 


Carriage  Speed 

Type 

magnetic  pick-up  to  counter 

(same  as  model  speed) 

Manufacturer 

Airpax 

Model 

Zero-Velocity  Pickup 

Wheel 

520-tooth 

Counter 

Hewlett  Packard 

Counter  Model 

5316B 

RPM 

Type 

magnetic  pick-up  to  counter 

Manufacturer 

Airpax 

Model 

Zero-Velocity  Pickup 

Wheel 

60-tooth 

Counter 

Hewlett  Packard 

Counter  Model 

5316B 

Drag  Dynamometer 

Manufacturer 

NSWCCD 

Model 

high  speed  500  lb 

Type 

full  bridge  strain  gauge  type 

Max  drag  force 

500  lb 

Max  lift  force 

500  lb 

Max  side  force 

500  lb 

Max  pitching  moment 

3200  ft-lbs 

Max  yawing  moment 

3200  ft-lbs 

Thrust  and  Torque  Dynamometer 

Manufacturer 

Kempf  &  Remmers 

Model 

R  58 

Type 

full  bridge  strain  gauge  type 

Rated  thrust 

+/-  1000  lb 

Rated  torque 

+/-  150  Ib-ft 

Max  rpm  in  water 

3600  rpm 

Overload  T&Q  allowed 

30% 

Signal  Conditioner 

Manufacturer 

Vishay 

(used  with  drag,  thrust,  and  torque  gage) 

Model 

2310 

Differential  DC  Amplifier 

Manufacturer 

Ectron 

(used  with  drag,  thrust,  torque,  rpm  and  speed)  Model 

751 ELN 
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Table  D.2.  Calibration  uncertainty  for  the  model  drag,  shaft  thrust,  and  shaft  torque  gages 
for  the  first  resistance  and  powering  test  series. 
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Note:  The  weights  used  for  these  calibrations  are  calibrated  every  2  years  to  a  tolerance  of  0.01% 


Table  D.3.  Calibration  uncertainty  for  the  model  drag,  shaft  thrust,  and  shaft  torque  gages 
for  the  second  resistance  and  powering  test  series. 
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Table  D.4.  Model  speed  measurement  uncertainty,  model  test  speed  range  ol 
for  the  first  resistance  and  powering  test  series. 


0.07%  0.07%  0.19%  0.10%  0.09%  avq  %  0.1 0% 


Table  D.5.  Model  speed  measurement  uncertainty,  model  test  speed  range  of  6- 1 8  knots 
for  the  second  resistance  and  powering  test  series. 
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Table  D.6.  Model  drag  measurement  uncertainty,  model  test  speed  =  6,10,16,18  knots, 
for  the  first  resistance  and  powering  test  series. 

10  knots 
16.878  ft/sec 


112  112 

133.487  133.423 

1.580  1.554 

1.18  1.16 


16  knots 

18  knots 

27.005  ft/sec 

30.3804  ft/sec 

RT  (lbs)** 

corr  RT 

RT  (lbs)** 

corr  RT 

(data  spot) 

(16  kn) 

(data  spot) 

(18  kn) 

#  of  spots 

average 

2*Stdev 

% 

25  25 

320.207  320.837 

3.768  3.121 

1.18  0.97 

26 

395.952 

4.410 

1.11 

26 

396.221 

4.515 

1.14 

Total  Average  % 

1.11 

**  each  data  spot  =  an  average  of  2000  samples  collected  over  a 
5  second  collection  time  at  a  rate  of  400  samples/sec. 


6  knots 


10.1268  ft/sec 

RT  (lbs)** 

corrected  RT 

(data  spot) 

(6kn) 

#  of  spots 

43 

43 

average 

47.482 

47.349 

2*Stdev 

0.595 

0.544 

%  1.25  1.15 
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Table  D.7.  Model  drag  measurement  uncertainty,  model  test  speed  =  6, 10,16, 18  knots, 
for  the  second  resistance  and  powering  test  series. 


6  knots 

10.1268  ft/sec 

10  knots 

16.878  ft/sec 

RT  (lbs)** 
(data  spot) 

corrected  RT 
(6kn) 

RT  (lbs)** 
(data  spot) 

corrected  RT 
(1  Okn) 

#  of  spots 

average 

2*Stdev 

% 

171 

47.688 

0.551 

1.16 

16  knots 
27.005  ft/sec 

171 

47.359 

0.521 

1.10 

38 

131.986 

1.300 

0.98 

18  knots 
30.3804  ft/sec 

38 

131.712 

1.148 

0.87 

RT  (lbs)** 
(data  spot) 

corrected  RT 
(1 6kn) 

RT  (lbs)** 
(data  spot) 

corrected  RT 
(18kn) 

#  of  spots 

1  4 

1 4 

62 

62 

average 

313.106 

312.790 

383.816 

383.846 

2*Stdev 

0.989 

0.922 

1.771 

2.094 

% 

0.32 

0.29 

0.5 

0.5 

Total  Average  %  0.70 


**  each  data  spot  =  an  average  of  2000  samples  collected  over  a 
5  second  collection  time  at  a  rate  of  400  samples/sec. 
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Table  D.8.  Model  thrust,  torque,  and  rpm  measurement  uncertainty,  model  test  speed=18  knots, 
for  the  first  resistance  and  powering  test  series. 


Table  D.9.  Model  thrust,  torque,  and  rpm  measurement  variation,  model  test  speed=18  knots, 
for  the  second  resistance  and  powering  test  series. 


75 


APPENDIX  E : 

UNCERTAINTIES  OF  FULL-SCALE  MEASUREMENTS 

The  relative  total  (precision  +  bias)  uncertainties  of  full-scale  measurements 
of  the  ship  speed,  propeller  rpm,  thrust,  torque,  and  shaft  horsepower  reported  for 
four  full-scale  trials  (USS  Boise  SSN764,  USS  Columbus  SSN762,  USS  Charlotte  SSN 
766,  USS  Memphis  SSN  691)  are  given  below 


Reported  uncertainties  of  full-scale  measurements 


ship 

speed 

rpm 

thrust 

torque 

SHP 

SSN  764 

0.6% 

1.9% 

4.1% 

0.3% 

1.9% 

SSN  762 

0.6% 

0.5% 

4.1% 

0.3% 

0.6% 

SSN  766 

0.3% 

0.2% 

2.2% 

0.9% 

0.8% 

SSN  691 

n/a 

0.2% 

2.2% 

1.4% 

0.9% 

Appreciable  variations  can  be  observed  in  the  foregoing  uncertainties.  Reasonable 
estimates  of  these  uncertainties  are  listed  below 


Uncertainties  of  full-scale  measurements  used  in  analysis 


speed 

rpm 

thrust 

torque 

SHP 

0.6% 

0.4% 

3.0% 

0.9% 

0.9% 

These  estimates  of  full-scale  measurement  uncertainties  are  used  in  the  present 
analysis. 
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